AD-A159  228  EXPLORATORV  DEVELOPMENT  ON  ADV. _ _ _ 

EXCITATION  AND  SIGN.  .  <U>  SOUTHHEST  RESEARCH  INST  SAN 
ANTONIO  TX  J  L  FISHER  ET  AL.  RPR  83  AFMAL-TR-83-4823 
UNCLASSIFIED  F33615-81-C-5836  F/G  14/2 


AD- 


A15&228 


AFWAL-TR-85-402  5 


EXPLORATORY  DEVELOPMENT  ON  ADVANCED  EDDY  CURRENT 
COIL  EXCITATION  AND  SIGNAL  PROCESSING 


Jay  L.  Fisher 
Dr.  R.  E.  Belssner 
James  E.  Doherty 


Southwest  Research  Institute 

6220  Culebra  Road 

San  Antonio,  Texas  7824 


April  1°A5 


Approved  for  public  release;  distribution  unlimited. 


MATERIALS  LABORATORY 

AIR  FORCE  WRIGHT  AERONAUTICAL  LABORATORIES 
AIR  FORCE  SYSTEMS  COMMAND 

VR I CHT- PATTERSON  AIR  PORCE  BASE,  OHIO  45433 


WJC  FILE  COPY 


s 


DTIC 

elect 

SEP  1  8 


T 


85  09  16  oOi 


mmcE 


when  Government  drawings,  specifications,  or  other  data  are  used  for  any  purpose 
other  than  in  connection  with  a  definitely  related  Government  procurement  operation, 
the  United  States  Government  thereby  incurs  no  responsibility  nor  any  obligation 
whatsoever/  and  the  fact  that  the  government  my  have  formulated,  furnished,  or  in 
any  way  supplied  the  said  drawings ,  specifications ,  or  other  data ,  is  not  to  be  re¬ 
garded  by  implication  or  otherwise  as  in  any  manner  licensing  the  holder  or  any 
other  person  or  corporation,  or  conveying  any  rights  or  permission  to  manufacture 
use,  or  sell  any  patented  invention  that  may  in  any  wav  he  mlated  thereto. 


This  report  has  been  reviewed  by  the  Office  of  Public  Affairs  (ASD/PA)  and  is 
releasable  to  the  national  Technical  Information  Service  (STIS).  At  ST IS,  it  will 
be  available  to  the  general  public,  including  foreign  nations. 


this  technical  report  has  been  reviewed  and  is  approved  far  publication. 


JAl^TTHOLLOWAY  / 

Nondestructive  Evaluation  branch 
totals  and  Ceramics  Division 


ram  rav  comma spot 

D.  :i.  FORNEY ,  JR.  t  Chl^ef^ 
Nondestructive  Evaluartun  Sfanch 
totals  and  Ceramics  Division 


mIf  your  address  has  changed,  if  you  wish  to  be  removed  from  our  mailing  list,  or 
it  the  addressee  is  no  longer  employed  by  your  organisation  p lease  notify  apuai  /mi  i  p 
I f-PAFB,  OS  45433  to  help  u a  maintain  a  current  smiling  list m. 


Copies  of  this  report  should  not  be  returned  unless  return  ia  required  by  security 
considerations,  contractual  obligations,  or  notice  on  a  specific  document. 


1*  RSPOMT  MCURITY  CLASSIFICATION 

UNCLASSIFIED 


JA  sscurity  classification  authority 


REPORT  DOCUMENTATION  PAGE 


1A.  RSSTRICTIVf  HARKINM 


OOCLASSIPICATIQM/OOWNORAOIMOSCMOOUL* 


A  Pf RPORMIMO  ORGANIZATION  IMPORT  NUMHIMI 


*A  NAM  OP  PS  RPORMIMO  ORGANIZATION 


AOOMN  tCtty.  SIMM  md  SIP  CaBM 

6220  Culebra  Road 
San  Antonio.  TX  78284 


MAMS  OP  PUNOINO/BPONSORINO 
ORGANIZATION 


t  0<STRISUTIOM/AVAI LABILITY  OP  IMPORT 

Approved  for  public  release;  distribution 
unlimited. 


S.  MOMITORIMO  ORGANIZATION  RIPORT  ttUMMRIII 

AR!AI.-TR-85-4025 


H  MAMS  OP  MOMITORIMO  ORGANIZATION 


tCUy.  1mm  md  ZIP  Cmdat 


Urlght-Patterson  AFB,  OH  45433-6533 


Sl  PROCURSMSNT  IMSTRUMSMT  IOSMTIPICATION  NUMBS  R 

F33615-81-C-5036 


10.  SOURCS  OP  PUNOINO  MOO. 


•  LS  MS  NT  MO. 


621 02F 


PROJOCT 

TASK 

WORK  UNIT 

NO. 

MO. 

NO. 

2418 

05 

16 

13.  PSRBOMAL  author*! 


13*  typb  op  roport 

ittna  in  — 

_ L 

13k.  TIMS  COVSRSO 

prom  9  Jon  81  TO 


1A  OATS 


•/ir-T.nnaf  Tvrrr— Him-: 

PI  i  I  1  H— — 


io  supplimsntary  notation 


COOAT  I  COOSS 


IS  «UMCT  TIRMI ICPIAMN  m  IWMR  S  hmmo  —PMM»Rft»y  MmO—n iMei 

Nondestructive  evaluation;  pulsed  eddy  current ‘..surface 
and  subsurface  flaws y  '  , 


ABSTRACT  iCmirw  mm  mwtrm  it  wwtn  «M  Mwlft  M  IMI  wtRl 

^  The  purpose  of  this  program  was  to  achieve  Improvements  In  eddy  current  performance 
by  analytical  and  experimental  exploratory  development  of  the  pulse  mode  of  coll 
excitation  as  a  means  of  evaluating  surface  and  near-surface  defects  In  low-conductivity 
materials  used  In  gas  turbine  engine  components. 

(  The  program  successfully  demonstrated  that  the  pulsed  eddy  current  approach  Is 
capable  of  detecting  flaws  as  small  as  0.010  Inch  by  0.005  Inch  In  a  T1-6-4  surface. 

In  addition  results  Indicate  that  this  approach  offers  a  promising  way  to  estimate 
the  depth  of  surface-breaking  flaws. 


30  OMTRtOUTlON/AVAIlAOlUTV  OP  ABSTRACT 
UNCLASSlPlSO/UNLIMlTtO  X  BAMS  AO  APT  □  OTIC  UOOI 


33a  NAMO  OP  RtSPONSlSLt  INDIVIDUAL 

Wn  O 

4.  A.  Hollows-  " 


31  ABSTRACT  SSCURITY  CLABBIPICATlON 

UNCLASSIFIED 


33k  rOLSPHONt  NWMBCR 


513)  255-5309 


S3*  OPPICO  SYMBOL 


PKEFACE 


This  is  cba  final  report  for  work  par  f  oread  bp  Southwest  Research  load  tut  a 
undar  AFWAL  Contract  Dunbar  F33615-81-C-S036,  Exploratory  Development  on 
Advanced  Eddy  Currant  Coil  Excitation  and  Signal  Processing.  Thia  project 
conalatad  of  an  investigation  of  natboda  to  laqprove  pulaad  addy  currant 
(PEC)  axaai nation  technology  for  flaw  detection  and  characterisation  in 
low  conductivity  naterlala  typically  uaed  in  gas  turbine  engine  compo¬ 
nents,  in  particular,  tltanlun  and  nickel-base  alloys. 


t  lion/ 

■  t  uty  Cadae 
»  aod/ar 

special 


IUU  or  C0NTEMT8 


!**• 


LIST  or  ILLUSTRATIONS .  1* 

!  1.  INTRODUCTION  AMD  StBMARY .  1 

1*1  Background . 1 

1.2  Objectives  . . 2 

1.3  Summary  of  Rooolto  .  .......  .  3 

1. A  Report  Sunnary . 3 

2.  BREADBOARD  SYSTEM  DEVELOPMENT  .  5 

2.1  Systam  Description  .....  .  5 

2.2  now  Specimens .  5 

2.3  now  Detection  Experiments .  7 

2.  A  Evaluation  Reeulte .  7 

2.3  Additional  Probe  Development  .  12 

3.  MODEL  DEVELOPMENT  AND  EXPERIMENTAL  VERIPICATION  .  19 

3.1  Introduction . 19 

3.2  Effects  of  Probe  Geometry  .  19 

3.3  Effects  of  Plow  Geometry .  23 

3.4  Comparison  with  Experiment-Zero  Crooning .  23 

3.3  Slot  Tip  Effects  -  Theory  . .  31 

3.6  Slot  Tip  Effects  -  Experiments  .  36 

3.7  Conclusion  ..............  .  39 

4.  ADVANCED  SYSTEM  DESIGN .  42 

4.1  Bread board  Development  .........  .  .  42 

4.2  Probe  Development  .  44 

4.3  Eve lust ion  Results  .....  .  ......  30 

3.  COMPARISON  WITH  C-W  RESULTS .  34 

6.  CONCLUSIONS  AND  RECOMSNDATIORS .  61 

7.  REFERENCES .  62 


APPENDICES 

A  Analytical  Models 

B  An  Alternative  Approech  for  Computing  Eddy  Currents! 

Case  of  the  Double-Layered  Plate 
C  Probe  Design  Studies 

D  Voltage  to  Current  Converter 

E  Flew  Specimens 


v 


LIST  OP  ILLUSTRATIONS 


A  PttlMd  Eddy  Current  Test  System  ........... 

Transmitted  Pulse  Shape  and  Frequency  Spec true  aa  a 
Function  of  Pulse  Rise  Tine  ...  . 

Detector  Scans  free  T1  6-4  Block  with  EDM  Notches  .  .  . 

Edfe  Resolution  Seans  froe  a  T1  6-A  Block  with  EDM 
Notches  .....  .  .  . 

Spatial  Resolution  Scans  froe  T1  6-4  Block  with  1.25  s 
.9  an  EDM  Notches  . 

Change  In  Inductance  of  a  Ferrite  Core  Coll  with 
Different  Shielding  Material  . 

Iaproved  Ferrite-Mound,  Not-to-Scale,  Ferrite-  and 
Copper-Shielded  Probe  with  Cylindrical  Syaaetry  .  .  .  . 

Toroidal-Wound  Probe  ....  . 

Coaparlsoo  of  Noraal  Component  of  Magnetic  Field  In  Air 

Simplified  Model  of  Pulsed  Eddy  Current  Probe  . 

Geometry  for  Detailed  Model  . 

Comparison  of  Calculated  and  Measured  Fields  froa  a 
Pulsed  Eddy  Current  Coll  In  Air  . 

Calculated  Signals  froa  Point  Flaws  at  Various  Depths 
Below  the  Surface  for  a  T  nsec  Gaussian  Pulse  . 

Point  Flaw  Signal  Peak  Amplitude  as  a  Function  of  Flaw 
Depth  for  Different  Pulse  Widths  ........... 

Confirmation  of  the  Existence  of  a  Liftoff-Insensitive 
Aaplltude  Croeslng  Point  . 

Calculated  Response  Function  for  a  .25  s  .15-aa 

(10  x  5-all)  Slot  . 

Calculated  Response  Function  for  a  1.0  x  .50-aa 

(AO  x  20-all)  Slot  . 


LIST  or  ILLUSTRATIONS  (Coat'd) 


Figaro  hit 

18  Calculated  and  Measured  Zero-Croesing  Times  aa  a 

Fuaetioo  of  Flaw  Depth  ....  .  31 

19  Flaw  Response  Baaed  on  aa  Integrated  Point  Flaw  Model  .  .  33 

20  Flaw  Response  Based  on  a  Crack  Tip  Model*  for  the  Sane 

Input  False  as  Deed  in  Figure  18  34 

21  Signal  Propagation  Paths  for  Open  end  Cloned  Cracks  ...  35 

22  Received  False  la  the  Absence  of  a  Flaw  (a)  Received 
Pulse  Vlth  Probe  Positioned  to  Maximize  the  Flew 

Signal  (b) .  37 

23  Flaw  Signal  Peak  Arrival  Time  Vs  Flew  Depth  aa  Measured 
(x'a)  Agrees  Vlth  the  wise  Required  for  ea  Eddy  Current 

False  to  Dlffoee  to  the  Creek  Tip  (solid  curve) .  38 

24  Ideal  Current  Waveforms  for  Achieving  Time -Domain 
Separation  of  Surface  Signals  and  Flaw  Depth 

Information .  40 

23  Modified  Pulsed  Bddy  Current  Test  System .  43 

26  Input  Voltage  Waveform  (Bottom  Trace)  and  Output  Current 

Waveform  (1  Ampere  per  Division)  of  the  Voltage-to- 
Current  Amplifier . . .  44 


27  Probe  Current  (Upper  Trace)  Follows  Input  Voltage 

(Lower  Trees)  Better  When  a  Probe  With  a  Higher 
Resonant  Frequency  ((b)*  14.3  MBs]  Probe  le  Deed 
than  When  a  Lower  Resonant  Frequency  Probe  ((e)* 


6  MBs}  is  used .  45 

28  The  Variation  of  Observed  Inductance  of  a  Single  Trans¬ 
mitter  Coll  aa  a  Function  of  Frequency  When  Measured 
with  Short  Leeds  and  vlth  2  Past  of  Different  Coaxial 

Cables .  47 

29  Actual  Current  Waveform  and  Power  Spectrum  Obtainable 

with  a  10-MHs  Probe  . . 48 

30  Magnitude  of  the  Fourier  Transform  of  the  Ideal  Trans¬ 

mitter  Current  Waveform  T  Is  the  Duration  of  the  Ramp 
Portion . 49 

31  No-Flaw  Recelvar  Signal  Used  In  a  Series  of  Flaw  Depth 

Measurement  Tests .  51 


vti 


LIST  or  ILLDSTIATIONS  (Coat'd) 


Deceiver  Plow  Ml  otto  Mo-Plow  Signal,  for  Dlf  farant  Oopth 
Motehos . . 

Invar  to  of  tho  Slop#  of  the  Decay  Portion  of  Plow 
Waveforms  Shown  la  Figure  32  . 


Conventional  Eddy  Currant  Sea no  of  Plows  A  Through  D 
on  tho  Surface  Plow  Detection  Standard  (E  la  Mot 
Detectable) . . . . 

Plan  View  lepreeentatloa  of  the  Vertical  (Hint nun 
Liftoff  Effect)  Couponont  of  tho  Inpedonce  Plane 
Display  for  the  Sane  Plows  as  in  Figure  3A  . 


Plan  View  iapreaantatlon  of  the  Vertical  (Mlnlnun  Liftoff 
Effect)  Component  of  the  Inpedonce  Plane  Display  for 
Plows  C,  H,  I,  and  J  . . . . 

Peak  Vertical  Signal  leaponee  free  Conventional  Eddy 
Current  Examination  ....  . 


Conventional  Eddy  Current  Scan  of  the  Sane  Plows  as  In 
Figure  36  Doing  an  Alr-learlng  Probe  Developed  by  SwKI 

for  MFC . . . 


Slnpllf led  Model  of  Pulsed  Eddy  Current  Probe 


Field  Calculated  0.023  an  la  Front  of  the  Model  Probe 
Shown  la  Figaro  39  ......  . 


Tine-Dependent  Plaid  Alongside  a  Probe  with  a  0.25-mn 
Thick  Shield  . 


Coapariooa  of  the  Analytic  Approxlaatloa  to  the  Exact 
Expression  of  Current  Density  . 


Geometry  Deed  in  Probe  Design  Studies 


Effect  of  the  Ferrite  Core  Site 


Effect  of  the  Coll-to-Shield  Distance 


Effect  of  the  Transmitter  Liftoff 


Effect  of  the  Copper  Shield 


■Vnr, 


m 


LIST  or  ILLUSTRATIONS  (Coat'd) 


Ul2£*  !y* 

48  Flaw  ittpoaM  Fttutlooi  for  Dlfforoatlol  Soaoora  ....  85 

49  Schoaatic  I  Hoot  ration  of  Tlaa-Oapaadant  Signal 

Faaturoa  ..........  .  87 

50  Sehoaotlc  of  Voltogo  to  Currant  Coavartar  .  90 


1.  INTRODUCTION  AMD  SIMtARY 


1.1  ltduromM 

Turbin*  engine  eoapoMot  ll(«^uuu|MMt  technology  now  taint  devel- 
op«d  under  a  "ret i rawer c-for-cauaa"  (ETC)  approach  requires  a  reliable 
net bod  (or  detecting  and  characterising  defects.  For  the  noat  part*  the 
present  high-resolution  techniques  for  defect  detection  applicable  to  rotat¬ 
ing  components  such  aa  discs  and  spacers  involve  single-frequency  eddy  cur¬ 
rent  examination.  As  evidenced  by  component  performance*  procedures  based 
on  these  techniques  have  effectively  screened  out  components  having  life- 
limiting  defects,  however*  these  simple  defect  detection  approaches  will 
not  satisfy  the  more  stringent  needs  of  RFC  programs.  New  methods  that 
provide  accurate  siting  of  defects  must  be  developed. 

PEC  examination  has  been  suggested  as  a  possible  approach  to  satisfy 
RFC  program  needs.  PRC  inspection  consists  of  generating  a  pulse  of  cur¬ 
rent  in  an  excitation  coil*  and  then  not  ins  the  characteristics  of  the 
signal  in  a  receiver  resulting  from  eddy  currents  excited  in  the  examina¬ 
tion  material.  This  is  in  contrast  to  conventional  eddy  current  examina¬ 
tion,  in  which  the  excitation  la  a  continuous  wave  of  one  or  more  discrete 
frequencies.  Another  difference  is  that  with  conventional  eddy  current 
examination,  aa  generally  practiced*  the  change  in  impedance  of  a  single 
(excitation)  coll  is  the  measured  parameter. 

There  are  several  possible  advantages  that  a  POC  eystmn  may  have  over 
a  conventional  eddy  current  system.  Plrst,  pulsed  operation  is  Inherently 
a  wideband  signal  excitation  technique*  and  thus  Includes  a  wide  rang*  of 
frequencies  without  the  need  to  explicitly  generate  each  frequency.  This 
fact  Is  of  Interest  because  electromagnetic  field  penetration  Into  a  con¬ 
ducting  material  varies  with  frequency  due  to  the  well-known  shin  effect* 
and  therefor*  nor*  Information  may  theoretically  be  obtaloed  by  using  sore 
frequencies.  Second,  there  Is  time-domain  information  available  in  PEC 
testing}  the  pulse  end  each  frequency  band  within  It  appear  at  specific 
Instances  In  time  after  traveling  through  the  examination  material*  because 
the  propagation  velocity  of  electromagnetic  fields  in  conducting  notarial 
Is  also  a  function  of  frequency.  Third*  pulsed  operation  coupled  with  a 
low  duty  factor  (puls*  repetition  rat*  times  pulse  duration)  for  the  trans¬ 
mitter  probe  allows  the  use  of  relatively  high  energy  excitation  pulses. 

For  example,  a  typical  duty  factor  of  0.1  percent  achieved  by  using 
1-microsecond  long  pulses  at  a  repetition  rat*  of  1000  pulses  per  second 
permits  high  energy  excitation  pulses  of  500  volts  naximun  and  S  amperes 
maximum  with  an  average  power  dissipation  of  less  than  2.)  watts. 

Although  the  technique  of  FtC  examination  has  existed  for  sosm  time, 
it  has  not  found  wide  application,  a  fact  that  Is  particularly  true  for  the 
low  conductivity  metals  of  interest  in  this  project,  historically,  only 
simple  approaches  have  been  used  to  extract  information  from  FEC  testing. 

Om  approach  ha*  been  to  simply  measure  the  peak  amplitude  from  a  differ¬ 
ential  receiver  coil  (1);  discontinuity  In  the  material  under  test  is 
indicated  when  the  ampTitude  varies.  Other  approaches  are  based  on  the 
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observation  that  thara  exist  certain  points  in  the  curve  of  output 
response  to  an  induced  pulse  that  do  not  change  significantly  when  one 
variable,  such  as  liftoff,  is  changed,  but  do  change  significantly  in  the 
presence  of  a  defect*  These  points  were  explained  by  Uaidellrh  in 
1956  0.  Other  suggested  sethods  have  involved  sampling  the  received 
signal  at  several  different  points  in  time.  The  aaplitudes  thus  obtained 
can  then  be  paraaetrlcally  combined  in  order  to  minimis#  the  effect  of 
varieties  of  features  not  of  interest  and  to  msximls*  the  signal  due  to 
features  of  interest  (for  example,  references  1,  4,  and  J).  This  approach 
is  motivated  by  a  desire  to  obtain  available  depth  sensitivity  similar  to 
ault if requeue y  eddy  current  testing,  and  takes  advantage  of  the  fact  that 
the  frequency  spectrum  varies  as  a  function  of  tiae  In  the  received  pulse. 
These  and  similar  approaches  that  involve  sampling  at  one  or  several  time 
or  magnitude  points  apparently  use  only  a  small  fraction  of  the  inherently 
wideband  signal  Information  that  ts  generated. 

An  additional  mot. vat  ion  for  new  research  into  PEC  testing  Is  that, 
despite  its  general  promise,  PEC  Instrumentation  has  generally  been  devel¬ 
oped  for  high  conductivity  and/or  low  magnetic  permeability  metals.  In 
such  material,  relatively  long  rlsetlae  (several  microsecond)  pulses  can 
be  used  because  the  propagation  time  of  magnetic  disturbance  Is  relatively 
long;  this  la  not  the  case  for  low  conductivity  nonmagnetic  materials  used 
In  gas  turbine  engine  disks.  Thus,  the  electronics  and  probes  for  these 
materials  wet  have  a  substantially  greater  bandwidth.  Solid  state  power 
electronics  has  improved  within  the  last  ten  years,  enough  so  that  precise 
control  of  pulses  In  this  frequency  range  at  reasonable  cost  Is  now  feasible. 

I .2  Objectives 

The  overall  objective  of  this  project  was  to  develop  *od  evaluate  a 
PEC  flaw  detection  capability  that  would  be  an  Improved  means  of  evaluating 
surface  and  near-surface  defects  in  low-conductivity  materials  used  In  gas 
turbine  engine  components.  The  specific  objective  wee  to  optimise  a  bread¬ 
board  pulsed  eddy  current  system,  including  the  probes,  instrumentation, 
and  signal  analysis  techniques.  A  specific  goal  was  the  reliable  detection 
of  o.ulo-inch  (0.2y-mn)  long  by  O.Wi-lnch  (O.I2>-m»)  wide  by  o.oui-lnch 
(U.WV-sm)  deep  sur face-breaking  flaws  that  are  of  interest  for  RPC  pur¬ 
poses.  An  ultimate  goal  Is  the  reliable  detection  and  siring  of  flaws  of 
this  site  and  larger,  for  both  surface-breaking  end  subsurface  flaws,  snd 
dlscrlelnatlon  of  such  defects  f roe  various  geometric  conf Igurstlons  such 
ss  el|M  snd  corners. 

The  princlpel  tasks  of  the  Svgl  approach  to  sccoepllsh  the  progrsn 
object  tees  were  as  follows: 

•k  1.  Analytical  Modeling  and  Probe  Optimisation 

Develop  a  simplified  model  for  predicting  PEC  response 
and  apply  to  probe  design  optimlsst ion.  Develop  s  more 
detailed  analytical  model  extending  classical  solutions 
for  coll  excitation  functions  to  the  transient  case,  for 
use  In  Identifying  signal  features  associated  with  oefect 
characteristics. 
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Task  II 


Assembly  of  an  Experimental  Pulsed  Eddy  Current  Breadboard 
Syatea 

Fabricate  reference  defect  specimens  for  system  evalua¬ 
tion.  Fabricate  optimized  probes,  assemble  a  pulse  gener¬ 
ator,  assemble  a  pulse  receiver,  and  develop  signal  pro¬ 
cessing  approaches.  Conduct  experiments  to  verify  the 
analytical  models.  Verify  the  performance  of  the  experi¬ 
mental  breadboard  system. 

Task  III.  Comparison  of  the  Pulsed  Breadboard  System  With  the  State- 
of-the-Art  Continuous-Wave  Eddy  Current  System 

Compare  the  results  of  the  pulsed  breadboard  system  with 
continuous-wave  instruments.  Provide  recommendations  for 
further  exploratory  development  of  the  PEC  approach. 

I  *  3  Summary  of  Results 

The  breadboard  PEC  system  successfully  detected  the  target-size 
flaws.  Initial  results  of  the  theory  correctly  predicted  the  field  distri¬ 
bution  from  sample  probes.  In  the  absence  of  a  flaw.  Therefore,  emphasis 
was  placed  on  the  possibility  of  using  unique  features  of  PEC  to  investi¬ 
gate  flaw  characterization.  It  was  found  that  the  simple  theoretical 
models,  obtained  by  integrating  point  flaw  model  fields  over  the  depth  of 
target-size  defects,  was  not  capable  of  quantitatively  predicting  the 
dependence  of  received  pulse  with  flaw  size.  Investigations  determined 
that  the  probable  cause  of  the  discrepancy  is  that  the  Integrated  point 
flaw  model  ignores  the  Kahn  effect  —  the  tendency  of  eddy  currents  to  be 
concentrated  at  the  tip  and  to  be  depleted  at  the  surface  of  a  surface- 
breaking  flaw.  It  was  realised  that  this  effect,  if  properly  detected  and 
used,  could  provide  information  about  flaw  depth,  since  the  time  that  it 
would  take  the  signal  from  the  tip-concentrated  eddy  currents  to  appear  at 
the  surface  would  be  characterlst lc  of  the  depth  of  the  source  of  concen¬ 
tration.  This  effect  was  modeled  by  the  use  of  subsurface  point  flaws. 
Experiments  confirmed  the  correlation  with  depth.  In  order  to  enhance  this 
type  of  measurement,  a  new  breadboard  system  was  constructed.  The  new  sys¬ 
tem  included  a  unique  pulser  that  allowed  precise  control  of  the  excitation 
coil  current.  The  pulse  was  arranged  to  minimize  eddy  current  excitation 
during  the  detection  period  of  the  pulse  cycle,  in  order  to  enhance  the 
signal-to-noise  ratio  of  the  depth  information.  The  existence  of  this  flaw 
depth  information  in  agreement  with  the  crack  tip  model  was  confirmed. 

1 •*  Report  Summary 

In  Section  2,  the  Initial  breadboard  system  is  described  and  success¬ 
ful  detection  of  target  flaws  is  shown.  Section  3  presents  the  model  devel¬ 
opment,  including  the  prediction  of  the  field  distribution  from  the  trans¬ 
mitter  probe  and  the  use  of  signal  characteristics  to  predict  flaw  depth. 
Comparison  with  experimental  data  is  presented.  It  is  shown  that  crack-tip 
eddy  current  concentration  (Kahn  effect)  must  be  included  in  the  model  in 
order  to  predict  flaw  depth  effects  and  to  agree  with  experimental  results. 
A  pulse  shape  to  enhance  the  depth  information  is  presented. 
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Section  4  describes  the  breadboard  system  ee  modified  to  toko  m»1h 
•dvuM«e  of  tbs  flow  depth  laf  ores  t  loo.  taper  lean  tel  results  show  agree* 
•eat  with  predictions  la  Section  3. 

Section  3  presents  tests  on  target  flews  using  conventional  continuous 
wave  eddy  current  examination  equipment.  It  la  shown  that  Its  detection 
capability  Is  comparable  to  that  of  the  Initial  breadboard  system,  but  Its 
use  in  defect  characterisation  Is  more  uncertain. 

Section  6  concludes  the  test  with  recoameodatlons  for  further  devel¬ 
opment  of  the  pulsed  eddy  current  method* 


2.  HEADBOARD  SYSTEM  DEVELOPMENT 


2.1  Eva  t—  Description 


Two  different  breadboard  lyttew  were  developed  durioc  thie  project. 
The  priaary  coaponents  of  each  systea  include  a  aechanical  scanner,  pulse 
generator,  signal  processor,  display  unit,  and  probes.  The  first  systea, 
described  in  this  section,  was  designed  to  test  beelc  concepts  in  probe 
design  and  flaw  detection.  Except  for  the  probes,  this  systea  was  con¬ 
structed  of  coaaercially  available  eleasnts.  The  scheaatic  layout  la  shown 
in  Figure  1.  The  Velooex  pulse  generator  la  capable  of  supplying  variable 
aaplltude  pulses,  with  peak  output  up  to  400  volts  into  s  noainal  200  oha 
load.  The  pulse  width  is  variable  froa  0.1  to  300  alcroeeconds. 

Signal  reception  froa  the  receiver  coil  is  accoapllshed  using  a  gen¬ 
eral  purpose  low-noise  differential  preaapllf ier,  with  a  linear  frequency 
response  froa  DC  to  50  HHs.  In  addition  to  aapllf lest ion,  the  preaapllfler 
any  be  used  to  provide  a  low-lapedance  (SO  oha)  load  to  the  receiver  probe. 
This  loed  any  be  thought  of  ea  providing  iapedance  notching.  In  the  tine 
does in  observations  that  have  been  node  in  this  project,  the  load  is  lapor- 
taot  in  reducing  the  circuit  <}  (quelity  factor)  of  the  probe-receiver  cir¬ 
cuit  to  osar  unity  so  that  its  response  to  a  rapid  voltage  change  is  well- 
daaped. 

A  significant  feature  of  the  systea  for  flaw  detection  la  the  boxcar 
averager  with  a  saapling  gate.  By  using  a  fixed-delay  gate  saaple  (a  fixed 
tlaa  interval  of  signal  duration  at  a  fixed  tlae  relative  to  the  pulsed 
trigger)  end  averaging  the  saaples  over  aany  pulse  repetitions,  unwanted 
noise  can  be  substantially  reduced.  The  integrating  function  of  the  lnstru- 
asnt  accumulates  DC  output  in  direct  proportion  to  the  average  of  the  sig¬ 
nal  within  the  gate  interval.  The  output  signal  is  thus  reduced  to  a  single 
value  of  enplitude  at  each  spatial  location.  The  gate  placement  can  be 
dictated  by  the  greatest  change  of  anplitude  on  the  waveforu  as  the  probe 
passes  over  the  artificial  flaws.  The  output  is  then  sent  to  an  E-Y 
recorder  or  to  a  single  channel  chart  recorder  for  archiving  and  further 
analysis. 

2.2  Flaw  Specleena 

Three  flaw  specinen  reference  standards  were  uanufactured  for  this 
project.  Detailed  data  sheets  in  Appendix  B  show  the  exact  flaw  layout 
with  the  specified  and  aeasured  flaw  sixes.  Each  standard  consists  of  s 
block  of  Tl-6-4  that  is  3  inches  (12.7  cn)  long  by  2  inches  (3.1  ce)  wide 
and  0.23  inch  (6.4  an)  thick.  Each  block  contains  a  series  of  EDM  notches. 
The  notches  are  all  0.003-inch  (0.13  ■)  wide.  They  range  in  length  froa 
0.030  inch  (1.23  an)  to  0.010  inch  (0.23  an)  and  in  depth  froa  0.030  inch 
(1.23  aa)  to  0.003  inch  (0.06  an).  The  defects  are  arranged  to  allow  six 
sets  of  continuous  datat 


(1) 

Equal  lengths,  incre 

aented  depths 

(2) 

Equal  depths,  iocrea 

anted  lengths 

»  v  r«j  -«  ' 


(3)  Different  lengths  end  depths,  with  eonstsnt  aspect  ratio. 

(4)  Parallel  pairs  of  notches  with  various  separation  distances. 

(5)  Notches  parallel  to  the  block  edge  at  different  distances  froa 
the  edge. 

(6)  Notches  perpendicular  to  the  block  edge  at  different  distances 
froa  the  edge. 

2.3  Flaw  Detection  Bxnerlaents 


The  probes  for  the  flaw  detection  experiaents  consisted  of  a  single 
transaltter  and  a  pair  of  aatched  receiver  coils.  Both  the  transaltter  and 
receiver  colls  were  encased  in  separate  copper  tubing  and  placed  coaxially 
next  to  each  other.  The  transaltter  and  receiver  colls  were  wound  on  sepa¬ 
rate  ferrite  cores  and  each  encased  In  s  6  aa  dlaaeter  copper  tubing  with 
a  0.9  aa  wall  thickness.  The  transaltting  coil  was  wound  using  SO  turns  of 
#32  gage  aagnet  wire. 

A  bipolar  pulse  froa  the  Velonex  pulser  was  used  to  excite  the  trans¬ 
altter.  The  risetiae  of  the  exaalnatlon  pulse  was  chosen  by  exaainlng  the 
frequency  content  of  the  excitation  current.  Figure  2  shows  the  effect  of 
pulse  shape  and  frequency  content  as  the  rise  tlae  of  the  current  pulse 
changed  froa  100  nanoseconds  to  20  nanoseconds.  The  excitation  voltage  of 
500  volts  was  applied  at  a  pulse  rate  of  300  cycles  per  second.  The  result¬ 
ant  pulse  width  ranged  froa  500  nanoseconds  to  1  aicrosecond.  A  risetiae 
of  ISO  nanoseconds  was  chosen  in  order  to  provide  a  flat,  wide  bandwidth 
signal.  This  pulse  has  a  reasonably  flat  frequency  response  out  to  1.2  MHz. 

The  receiver  colls  consisted  of  a  pair  of  ISO  turn,  #40  coils  wound 
on  a  single  ferrite  core  to  fora  a  differential  receiver  probe.  The 
receiver  colls,  aa  well  as  the  transaltting  coll,  were  encased  in  the  6  aa 
copper  tubing. 

2.4  Evaluation  Results 

The  probe  described  above  was  scanned  across  the  surface  of  the  aanu- 
factured  flaw  speclaens*  Figure  3  indicates  test  results  obtained  using  a 
tltanlua  6-4  surface  flaw  speciaen.  The  sasllest  flaw  detectable  is  2S0 
alcrons  long  by  12S  aicrona  deep.  In  each  case,  the  probe  was  oriented 
such  that  a  line  between  the  transaltter  and  receiver  colls  was  parallel 
to  the  notch  length. 

Figure  4  Indicates  the  edge  resolution  capability  of  the  PBC  probe, 
using  data  froa  the  edge  resolution  standard.  This  standard  has  1.25  an 
long  by  0.9  an  deep  notches  oriented  perpendicular  to  the  edge  of  the  plate 
at  various  distances  froa  the  edge  (as  described  in  Appendix  E).  The  clos¬ 
est  distance  at  which  a  notch  was  detected  was  1.25  aa  froa  the  edge  of  the 
plate. 


■.V.v’VWV 


(Cont’d).  Transaicced  Pulse  Shape  and  Frequency  Spectrua  aa  a  Function  of  Pulse  Rise 
diaaeter  probe  with  a  0.9-aa  wall  was  used.  For  low  rise~tiae  pulses,  the  probe 
resonant  peak. 


Figure  4.  Edge  Resolution  Scans  froa  a  Ti  6-4  Block  with 
EBM  Notches.  Notches  are  1250  pa  long  and  250  ua  deep. 
These  data  are  taken  using  the  saae  probe  used  in  Figure 
All  notches  ere  detected. 
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Figure  5  indicates  the  results  of  the  spatial  resolution  test 
obtained  with  the  saae  PEC  probe.  The  specimen  scanned  contains  sets  of 
EDM  notches,  1.25  aa  long  by  0.9  aa  deep  by  0.125  aa  wide,  arranged  paral¬ 
lel  to  each  other  at  various  apacings,  aa  indicated  in  the  figure.  Aa  can 
be  seen,  it  is  not  possible  to  separate  two  flaws  that  are  closer  than 
2.5  sai  with  this  probe.  Probes  with  saaller  disasters,  however,  can  sepa¬ 
rate  flaws  spaced  as  close  as  1.25  aa.  These  are  described  in  the  next 
section. 

2.5  Additional  Probe  Developaent 

Early  probes  were  designed  with  siaple  geoaetries  to  aatch  the  aodel- 
lng  capability  described  in  the  following  section.  These  probes  demon¬ 
strated  that  the  basic  flaw  detection  goals  of  the  prograa  could  be  aet. 
However,  the  spatial  resolution  capability  of  these  probes  was  only  margi¬ 
nal  .  Therefore,  additional  probe  developaent  was  concerned  with  laproving 
this  capability.  The  basic  approach  was  to  minimize  the  electroaagnetlc 
footprint  of  the  transalttlng  coll  on  the  inspection  surface.  This  aay  be 
accomplished  by  asking  the  probe  saaller,  of  course,  and  by  using  a  probe 
design  with  better  focusing  capability.  Both  of  these  approaches,  espe¬ 
cially  the  former,  tend  to  reduce  the  magnetic  field  aaplitude  in  the 
desired  Inspection  area. 

One  way  to  provide  a  saaller  footprint  without  sacrificing  size  is 
to  improve  shielding.  Figure  6  shows  the  effect  of  different  shielding 
material  and  lack  of  shielding.  The  probe  consists  of  50  turns  of  #36 
gage  wire  wound  around  the  ferrite  core.  The  apparent  Inductance  aa  a 
function  of  frequency  is  shown  for  the  coil  unshielded,  and  shielded  with 
copper,  alualnua,  steel,  and  ferrite.  The  inductance  produced  is  directly 
proportional  to  the  arount  of  energy  created  in  the  magnetic  field  gener¬ 
ated  by  the  probe.  Since  the  magnetic  energy  density  is  proportional  to 
the  square  of  the  aagnetlc  field  intensity,  a  reduction  in  magnetic  field 
volusw  and  proportional  Increase  in  flux  density  causes  an  Increase  in 
total  energy  and  therefore  Inductance.  A  decrease  in  Inductance  aay  be 
caused  by  an  increased  magnetic  volume,  but  also  by  a  reduction  of  inten¬ 
sity  due  to  increased  generation  of  eddy  currents  in  the  probe.  Copper, 
steel,  and  aluadnua  provide  shielding  by  generating  eddy  currents  that  can¬ 
cel  the  primary  field.  Since  ferrite  Is  s  poor  conductor,  it  does  not  gen¬ 
erate  an  opposing  aagnetlc  field;  rather  it  provides  a  low  reluctance  path 
for  the  aagnetlc  field  lines  exterior  to  the  probe  core  between  the  top  of 
the  core  and  the  inspection  material.  The  figure  shows  that  at  the  rela¬ 
tively  high  testing  frequencies  of  interest  in  this  project,  eddy  current 
losses  in  copper,  alualnum,  and  steel  cause  a  substantial  loss  in  magnetic 
field  energy,  even  though  the  aagnetlc  volume  must  be  severely  reduced. 
These  losses  Increase  as  the  frequency  increases.  On  the  other  hand,  the 
ferrite  shielding  increases  the  Inductance  by  a  factor  of  approxlaately  1.8 
over  the  entire  frequency  range  of  interest.  It  is  concluded  that  ferrite 
shielding  is  quite  useful  in  concentrating  the  sugnetic  field  energy,  and 
that  highly  conducting  aaterlal,  while  of  use  for  low  frequency  eddy  cur¬ 
rent  testing  is  not  as  desirable  at  the  higher  frequencies  appropriate  to 
examination  of  poorly  conducting  alloys. 


v.v:-v. : 
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In  order  to  obtain  quantitative  Measurements  oi  the  focusing  ability 
of  copper  veraua  ferrite  shielding,  Measurements  were  taken  of  the  normal 
magnetic  field  produced  after  passage  through  two  different  thicknesses  of 
titanlua  for  two  probes  with  the  saae  else  ferrite  core.  The  amplitude  of 
the  normal  component  of  the  field  was  Measured  by  a  smsII  (1  an  dlaMeter) 
sensing  coll,  as  a  function  of  distance  along  one  axis  at  the  surface  oppo¬ 
site  the  transMitter,  for  thickness  of  1.3  and  3.2  mm.  Results  for  the 
footprint  sise,  ss  indicated  by  the  distance  over  which  the  SMplltude  Is 
at  least  half  the  uaxluuu  amplitude,  are  shown  in  Table  1.  It  is  seen 
that  the  ferrite  shielding,  in  addition  to  providing  greater  field  strength 
has  a  footprint  approxiMately  the  sane  as  the  copper-shielded  probe  at  the 
eMail  thickness,  but  which  Increases  only  29  percent  coapared  to  47  percent 
through  the  entire  thickness  of  titaniuM,  so  that  the  final  footprint  is 
20  percent  less  than  that  of  a  eluilarly  sited  copper-shielded  probe. 


Table  1 

FULL  WIDTH  OF  HALF  MAXIMUM  NORMAL  FIELD  IN  THROUGH  TRANSMISSION 


Ti  Thickness  Copper-Shielded  Probe 


Ferrite-Shielded  Probe 


0.030-  (1.27  mm)  3.0  mm 

0.123-  (3.18  mm)  4.4  mm 


2.6  ms 


3.6  MM 


Two  additional  iaprovewnts  have  been  Made  to  the  ferrite  shielded 
probe.  It  has  been  found  that  a  layer  of  copper  shielding  placed  outside 
of  the  ferrite  can  be  used  to  reMove  stray  fields  reMSinlng,  and  thus 
reduce  coupling  to  the  receiver  probe.  For  a  sufficiently  thick  ferrite 
shield,  thts  additional  layer  has  negligible  effect  on  the  probe  iape- 
dance,  but  it  does  reduce  the  coupling  to  an  unshielded  receiver  probe. 

The  aMount  of  reduction  depends  on  receiver  placement. 

A  second  lMproveMent  Is  the  placeMent  of  a  copper  disk  around  the 
ferrite  core  between  the  core  and  the  shielding  at  the  exaMlnatlon  surface 
side  (Figure  7  shows  the  coMpleted  probe  design).  This  disk  is  used  to 
increase  the  Magnetic  reluctance  of  the  air  path  between  the  shielding  and 
the  core,  in  order  to  MinlMite  the  nuMber  of  leakage  field  lines  that  pass 
froM  the  core  to  the  shield,  tangential  to  the  inspection  surface,  without 
penetrating  the  surface.  In  order  to  prevent  eddy  currents  in  the  copper 
disk  froM  also  opposing  the  normally  oriented  Magnetic  field,  a  vertical 
slit  is  cut  through  the  entire  height  of  the  disk  at  one  circuMferential 
location. 


A  possible  problea  with  the  use  of  the  spilt  ring  conducting  disk  is 
that  a  small  gap  aay  have  a  large  enough  capacitance  to  allow  significant 
eddy  current  density  to  develop  in  opposition  to  the  primary  field.  This 


problsm  U  |r««t«r  «t  higher  (r«qwscl«i,  line*  the  capacitive  reactance 
U  Inversely  preport lonel  to  frequency .  However,  neasureneots  taken  with 
this  probe  design  show  an  increase  in  the  oornal  component  of  the  nagnetic 
field  of  1U  percent  through  a  3.2-aa  thick  titanium  plate,  conparnd  to  a 
sintlar  probe  without  the  copper  dish. 

Another  type  of  transmitting  coil  that  was  evaluated  used  a  toroidal 
gap  core.  The  transmitter  geometry  la  illustrated  in  Figure  8.  The 
toroidal  cores  had  0. 325-inch  outside  diameter,  0. 187-inch  Inside  diameter, 
and  were  0.125-loch  thick.  Cap  distances  ranged  from  0.050  to  0.100  inch. 
Figure  9  shows  the  normal  component  of  magnetic  field  for  a  coaxial  probe 
encased  in  6-nm  diameter  copper  tubing,  and  the  same  component  across  the 
air  gap  of  a  toroidal  gap  core.  It  was  concluded  that  the  shielded  coaxial 
probe  offers  better  focusing  than  the  toroidal  gap  coll,  although  with  cop¬ 
per  shielding,  the  field  intensity  of  the  coaxial  probe  is  approximately 
one-half  that  of  the  toroidal  gap  coll.  With  the  use  of  ferrite  shielding, 
this  disadvantage  Is  ameliorated.  Therefore,  It  was  decided  not  to  pursue 
this  design. 
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Figure  7.  Iaproved  Ferrite-Wound ,  Not-to-Scale,  Ferrite 
and  Copper-Shielded  Probe  with  Cylindrical  Syeaetry. 
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Figure  8*  Toroidal-Wound  Probe.  This  probe  showed 
less  focusing  than  cylindrical  probes. 


Figure  9.  Comparison  of  Normal  Component  of  Magnetic  Field  in  Air.  The 
top  field  profile  is  obtained  from  a  shielded  coaxial  probe,  while  the 
bottom  is  from  a  toroidal  gap  coil.  As  shown  by  the  insets,  similar 
current  pulses  were  used  for  both  coils.  The  shielded  coaxial  probe 
exhibits  less  spatial  broadening. 


3.  MODEL  DEVELOPMENT  AMD  EXPERIMENTAL  VERIFICATION 


3.1  Introduction 


Th«  development  of  a  Mathematical  Model  of  a  pulsed  eddy  current 
aystem  was  approached  In  two  steps*  First,  the  simplified  model  Illus¬ 
trated  In  Figure  10  was  developed  to  explore  the  effects  of  probe  geometry 
on  the  fields  produced  In  the  absence  of  a  conducting  test  piece.  The 
second  step  Included  a  generalisation  of  the  simple  model  to  Include: 

(1)  A  finite  coll  winding  thickness, 

(2)  The  approximate  effect  of  a  ferrite  coll, 

(3)  A  separate,  differential  receiver  coll, 

(4)  The  change  in  receiver  response  caused  by  the  presence  of  an 
unflawed  apeciman,  and 

(3)  Additional  changes  caused  by  the  presence  of  a  point  flaw. 

Parameters  included  In  the  detailed  model  are  those  shown  in  Fig¬ 
ure  11.  A  more  complete  deecription  of  both  models,  along  with  the 
results  of  some  exploratory  calculations,  is  given  In  Appendix  A.  An 
alternate  model,  which  required  too  much  computer  time  for  pulsed  eddy 
current  analyses  but  which  may  be  useful  for  other  purposes,  is  described 
In  Appendix  B. 

3.2  Effects  of  Probe  Geometry 

As  e  test  of  the  model  of  the  transmitter  coll  and  shield,  calcula¬ 
tions  of  the  magnetic  field  In  air  were  performed  and  compared  with  experl' 
mental  data  obtained  by  scanning  a  small  Hall-effect  sensor  radially 
across  the  coll.  The  result,  which  is  shown  in  Figure  12,  verifies  the 
validity  of  prediction  of  fields  produced  by  a  shielded  excitation  coil. 

The  detailed  model  was  then  applied  to  parametric  studies  of  the 
effects  of  probe  geometry  on  the  eddy  current  distribution  In  a  conducting 
material  and  on  probe  response  to  a  hypothetical  point  flaw.  The  results 
of  these  calculations,  which  arc  summarised  in  Appendix  C,  show  that  the 
effects  of  probe  geometry  are  similar  to  those  observed  with  continuous 
wave  eddy  current  systema.  Thus,  for  example,  the  calculations  show  that 
the  eddy  current  distribution  is  spatially  broadened  and  increaaea  in 
intensity  as  the  transmitter  coil  diameter  is  increased.  Small  changes  in 
the  coll-to-shield  distance  produce  minor  variations  In  the  distribution, 
though  removal  of  the  shield  leads  to  a  more  diffuse  pattern  of  current 
flow  in  the  material.  Use  of  a  differential  sensor  leads  to  improved  flaw 
resolution,  as  expected,  with  the  maximum  signal  strength  occurring  when 
the  flaw  is  directly  below  the  transmitter  coil  winding.  In  summary, 
then,  the  model  reveals  that  the  influence  of  probe  geometry  on  pulsed 
eddy  current  flaw  response  is  essentially  the  same  as  in  the  continuous 
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Figure  10.  Simplified  Model  of  Pulsed  Eddy  Current  Probe.  This 
model  was  used  to  evaluate  the  field  in  front  of  and  alongside 
the  probe.  The  finite  current  coll  is  modeled  as  a  current  sheet 
of  infinitesimal  thickness. 


wave  case.  It  follows  that,  as  far  as  the  geometrical  aspects  of  probe 
design  are  concerned,  experience  with  continuous  wave  probes  is  directly 
applicable  to  the  design  of  pulsed  eddy  current  probes. 


The  distinguishing  feature  of  pulsed  eddy  current  flaw  response  is, 
of  course,  the  time  dependence  of  the  signal  which,  hopefully,  contains 
information  pertaining  to  flaw  geometry.  Thus,  while  calculations  of  the 
influence  of  probe  geometry  were  needed  to  determine  possible  deviations 
from  the  continuous  wave  case,  the  principal  reason  for  undertaking  pulsed 
eddy  current  model  development  was  to  explore  features  in  time-dependent 
flaw  signals  and  their  relationship  to  flaw  geometry. 

3.3  Effects  of  Flaw  Geometry 

The  results  of  some  early  calculations  of  flaw  signals  are  shown  in 
Figures  13  and  14.  Figure  13  shows  calculated  signals  from  point  flaws  at 
various  depths  below  the  surface  plotted  as  a  function  of  time  (in  usee) 
and  normalized  to  unity  at  peak  value.  The  important  point  illustrated 
here  is  that  flaw  signal  peak  arrival  times  increase  with  flaw  depth  and 
therefore  provide  a  measure  of  depth.  However,  for  point  flaw  depths 
corresponding  to  crack  tip  depths  of  interest,  peak  arrival  times  are  of 
the  order  of  tens  of  nanoseconds.  Thus,  if  depth  information  is  to  be 
obtained  from  peak  arrival  time  measurements,  input  pulse  widths  or  pulse 
decay  times  must  be  extremely  short,  probably  of  the  order  of  a  few 
nanoseconds . 

Another  characteristic  of  point  flaw  signals  is  illustrated  in 
Figure  14.  This  is  a  plot  of  flaw  signal  peak  amplitude  as  a  function  of 
flaw  depth  for  three  different  pulse  widths.  The  main  point  to  note  here 
is  that  predicted  signal  amplitudes  decrease  by  about  a  factor  of  100  in 
the  first  0.1  mm  of  depth  below  the  surface.  Thus,  based  on  these  two 
sets  of  point  flaw  calculations,  one  might  expect  that  signals  of  interest 
are  very  weak  compared  to  signals  from  near-surface  reflections  and  that 
the  arrival  times  of  interest  are  of  the  order  of  tens  of  nanoseconds. 

However,  experimental  results  for  surface-breaking  slots  of  various 
depths,  which  will  be  presented  in  the  next  section,  show  a  very  definite 
depth  dependence  even  though  the  input  current  waveform  was  much  broader 
than  that  used  in  the  calculation  shown  in  Figure  13.  This  observation, 
which  seemed  to  be  in  conflict  with  model  calculations,  called  for  more 
detailed  comparisons  between  theory  and  experiment  and  led,  eventually,  to 
a  reexamination  of  one  of  the  basic  assumptions  inherent  in  the  model. 

3.4  Comparison  With  Experiment-Zero  Crossing 

A  series  of  experiments  was  carried  out  in  order  to  compare  measured 
flaw  response  with  theoretical  predictions.  The  parameter  measured  was 
zero-crossing  time,  that  is,  the  time  relative  to  an  arbitrary  beginning 
of  pulse  time  that  the  received  signal  first  crosses  below  zero  amplitude. 
This  measurement  was  chosen  for  two  reasons.  First,  it  provides  an  indica 
tion  of  the  peak  signal  arrival  time,  but  is  more  precisely  measured  than 
the  peak  time  itself  because  the  slope  of  the  signal  is  in  general  greater 
at  the  zero-crossing  point.  Second,  previous  theoretical  and  experimental 
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Figure  13.  Calculated  Signals  from  Point  Flaws  at  Various 
Depths  Below  the  Surface,  for  a  1  nsec  Gaussian  Pulse.  Fo 
corresponding  crack  tip  depths  of  interest,  peak  signal 
arrival  time  is  of  the  order  of  10  nsec. 
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Figure  14.  Point  Flaw  Signal  Peak  Amplitude  as  a  Function  of  Flaw 
Depth,  for  Different  Pulse  Widths.  Predicted  signal  amplitude 
decreases  by  a  factor  of  approximately  100  in  the  first  .1  mm  of 
depth,  indicating  very  weak  signals  compared  to  near-surface 
reflections  for  this  model. 
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results  have  shown  that  there  Is  in  general  an  amplitude,  the  time  of 
which  is  sensitive  to  the  presence  of  flaws  but  not  to  variation  in  lift¬ 
off  (2).  This  second  factor  can  be  predicted  by  assuming  that  liftoff 
variation  is  equivalent  to  a  change  in  receiver  probe  resistance  and  induc¬ 
tance.  The  simple  theory  does  not  guarantee  that  the  amplitude  predicted 
will  actually  be  in  the  range  of  amplitude  available.  However,  it  was 
experimentally  confirmed  that  this  point  exists  for  the  probes  in  the  PGC 
breadboard  system  on  a  titanium  sample,  as  shown  in  Figure  15.  It  is 
observed  from  this  figure  that  this  amplitude  crossing  point  is  close  to 
the  zero-crossing  point.  With  this  justification,  a  series  of  experiments 
measuring  zero-crossing  time  was  carried  out.  A  copper-shielded  trans¬ 
mitter  probe,  with  differential,  shielded  receiver  coils  was  used.  The 
results  will  be  presented  after  a  discussion  of  the  theoretical 
calculations. 

In  Figure  15(a),  the  receiver  probe  signal  over  Ti-6-4  sample  with 
no  flaw,  for  three  different  values  of  liftoff,  is  shown.  Note  that  at 
0.5  amplitude  (vertical)  divisions,  the  crossing  time  is  approximately 
independent  of  liftoff.  In  Figure  15(b),  the  receiver  signal  over  a 
1.8  mm  deep  flaw  is  shown,  again  for  three  different  values  of  liftoff. 

The  liftoff-insensitive  crossing  point  exists  and  is  at  the  same  near-zero 
amplitude  but  occurs  160  nanoseconds  earlier. 

To  make  direct  comparisons  between  theory  and  the  zero-crossing  data 
described  above,  it  was  necessary  to  modify  the  analytical  model  to 
account  for  the  geometry  of  the  flaws  used  in  the  experiment.  In  making 
this  modification,  it  was  necessary  to  introduce  additional  approximations 
because  the  basic  theory  makes  use  of  the  assumption  that  the  electric 
field  inside  a  flaw  is  proportional  to  the  unperturbed  field  at  the  center 
of  the  flaw.  This  can  be  rigorously  justified  only  in  the  case  of  small 
subsurface  flaws  such  that  the  unperturbed  field  does  not  vary  signifi¬ 
cantly  over  the  surface  of  the  flaw.  However,  for  an  extended  flaw  in  a 
nonuniform  field,  one  might  expect  the  interior  field  to  vary  from  point 
to  point  in  a  way  that  is  approximately  proportional  to  the  spatial  varia¬ 
tion  of  the  unperturbed  field.  Thus,  as  a  first  approximation  to  a  theory 
for  an  extended  flaw,  the  point  flaw  solution  can  be  simply  integrated 
over  the  volume  of  the  flaw.  The  resulting  flaw  response  function  that 
depends  on  flaw  dimensions  is  multiplied.by  a  scattering  factor  a,  which 
is  a  scalar  approximation  to  the  matrix  a  defined  in  Appendix  A.  Because 
there  is  no  theory  that  allows  one  to  calculate  the  value  of  a  for  a 
rectangular  slot,  it  was  treated  as  an  adjustable  parameter  in  making 
comparisons  with  the  experiment. 

However,  in  experimental  data,  the  flaw  response  function  alone  is 
not  observed,  but  rather  the  sum  of  the  signal  with  no  flaw  and  the  change 
caused  by  a  flaw  are  observed.  Thus,  to  compare  with  the  experiment,  one 
must  add  a  no-flaw  component  to  the  calculated  signal.  Because  the  no-flaw 
signal  is  liftoff  dependent,  this  amounts  to  incorporating  liftoff  depend¬ 
ence  as  part  of  the  model. 
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Figure  15.  Confirmation  of  the  Existence  of  a  Liftoff- 
Insensitive  Amplitude  Crossing  Point.  At  top  is  the 
no-flaw  signal  for  3  different  values  of  liftoff.  At 
the  bottom  is  the  flaw  signal  for  3  different  values  of 
liftoff,  with  the  crossing  point  occurring  .8  units 
earlier. 


A  mathematical  model  of  the  response  of  a  pulsed  eddy  current  probe 
to  liftoff  variation  was  therefore  developed.  The  theory  is  based  on  a 
time-dependent  generalisation  of  Auld's  reciprocity  theorem  (6);  the  cal¬ 
culation  proceeds  In  much  the  same  way  as  in  the  continuous  wave  eddy 
current  case.  The  principal  result  is  the  following  expression  for  the 
change  in  the  Impulse  response  function  as  the  probe  is  moved  from  an 
Infinite  liftoff  distance  to  liftoff  distance  t. 
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where  Rj  and  R£  are  the  radii  of  excitation  and  differential  receiver 
coils,  respectively;  Jg  and  are  Bessel  functions  of  the  first  kind  of 
order  0  and  1,  respectively;  and  t ,  the  scaled  time,  is  given  by 
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The  response  of  a  pulsed  eddy  current  probe  is  the  convolution  of 
the  excitation  pulse  with  an  Impulse  response  function  which  is  the  sum  of 
r0,  the  response  in  air  in  the  absence  of  a  specimen,  fiTg,  the  change 
caused  by  an  unflawed  specimen  at  liftoff  £,  and  5Tf,  the  change  caused  by 
a  flaw.  As  was  noted  earlier,  the  magnitude  of  6Tf  relative  to  the  other 
components  of  the  Impulse  response  function  Is  unknown  because  the  calcula¬ 
tion  of  is  based  on  a  small  flaw  theory  and  contains  an  adjustable 
parameter  to  approximately  account  for  finite  flaw  dimensions.  Therefore, 
the  change  in  response  caused  by  the  specimen  and  flaw  can  be  described  by 


sr  -  «rt  +  a6rf 


where  a  is  an  adjustable  multiplier. 

In  the  initial  set  of  liftoff-dependent  calculations,  an  exponen¬ 
tially  damped  sine  wave  was  used  as  the  excitation  current,  and  a  was 
varied  in  an  attempt  to  obtain  qualitative  agreement  with  experimental 
data.  Typical  results  are  shown  in  Figures  16  through  17.  The  dotted 
curves,  which  *>re  the  same  in  both  figures,  are  predictions  of  probe 
response  in  the  absence  of  flaws.  The  dashed  curves  are  flaw-induced 
changes,  and  the  solid  curves  are  the  sums  of  the  dotted  and  dashed 
curves.  The  parameter  DT  is  one-half  the  period  of  the  excitation  pulse 
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Figure  16.  Calculated  Response  Function  for  s 
.25  x  .15-aa  (10  x  5-ail)  Slot. 


Figure  17.  Calculated  Response  Function  for  a 
1.0  x  .50-an  (40  x  20-ail)  Slot. 


in  microseconds,  DEPTH  Is  the  depth  In  ails  of  a  rectangular  surface  slot 
with  length  •  2  x  DEPTH,  and  FACTOR  Is  a  nuaber  proportional  to  the  adjust 
able  parameter  a. 

The  effect  such  calculations  were  designed  to  explain  is  the  experl- 
aentally  observed  shift  in  the  first  zero  crossing  as  a  function  of  flaw 
size.  Fro*  the  figures.  It  can  be  seen  that  the  theory  attributes  this 
shift  to  an  Increase  in  the  aagnltude  of  the  flaw  signal  relative  to  the 
no-flaw  signal  as  flaw  size  Increases.  Thus,  for  exaaple,  in  Figure  16, 
the  flaw  signal  (dashed  curve)  is  saall  compared  to  the  no-flaw  signal 
(dotted  curve);  and  the  zero-crossing  shift,  the  difference  between  the 
dotted  and  solid  curves,  is  saall.  In  Figure  17,  on  the  other  hand,  the 
flaw  signal  is  auch  larger  coapared  to  the  no-flaw  signal;  and  the  zero- 
crossing  shift  is  correspondingly  greater.  Significantly,  these  data  indl 
cate  that  the  zero  crossing  is  alaost  entirely  deterained  by  the  relative 
aagnltudes  of  the  flaw  and  no-flaw  signals;  there  is  no  significant  change 
In  the  shape  of  the  flaw  signal  as  flaw  size  Increases. 

Calculated  zero-crossing  data  are  plotted  in  Figure  10  along  with 
experlaental  data.  The  calculated  tlaes  are  shifted  approximately 
0.3  usee  to  agree  with  the  aeasured  zero-crossing  tlae,  about  1.3  i-sec, 
with  no  flaw  present.  As  can  be  seen  froa  this  figure,  both  theory  and 
experlaent  show  an  increase  in  zero  crossing  tlae  with  decreasing  flaw 
depth,  though  the  functional  foras  of  predicted  and  aeasured  relationships 
differ  significantly. 

To  deteralne  If  this  difference  is  caused  by  a  poor  choice  of  the 
adjustable  paraaeter  a  that  fixes  the  flaw  response  aaplitude,  we  allowed 
a  to  vary  froa  about  one-half  the  value  used  here  to  a  aaxlaua  of  four 
tlaes  this  value.  In  all  cases,  the  predicted  curve  was  concave  upward, 
as  in  Figure  18,  In  consistent  dlsagreeaent  with  the  concave  downward 
curve  deterained  froa  experlaent. 

One  possible  reason  for  this  discrepancy  is  that  the  model  does  not 
Include  the  probe  response  in  air  in  the  absence  of  a  speciaen,  l.e.,  the 
tera  To  mentioned  above.  If  this  tens  were  added  to  the  calculation,  the 
dotted  curves  in  Figures  16  and  17  would  change  shape,  which  could  result 
in  a  different  behavior  for  the  zero  crossing  tlae.  However,  the  fact 
that  slailar,  concave  upward  curves  like  that  in  Figure  18  have  been 
obtained  for  different  excitation  waveforms,  which  also  changes  the  shape 
of  the  dotted  curve,  would  seea  to  Indicate  that  the  problea  lies  with  the 
flaw  signal  rather  than  the  no-flaw  response. 

In  particular,  the  fact  that  this  aodel  predicts  no  significant 
change  in  the  shape  of  the  flaw  signal  as  flaw  size  is  varied  could  very 
well  be  an  artifact  of  the  saall  flaw  approxlaat ion  which  foras  the  basis 
for  this  aodel  of  a  rectangular  slot.  Froa  other  calculations  of  the 
behavior  of  continuous-wave  eddy  currents  at  corners  and  around  crack 
tips,  one  expects  that  there  is  a  depletion  of  current  density  near  the 
slot  on  the  surface  of  the  speciaen,  and  an  enhancement  of  current  density 
at  the  slot  tip.  Neither  of  these  so-called  "Kahn  effects"  are  predicted 
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Figure  18.  Calculated  and  Measured  Zero-Crossing  lines  as  a 
Function  of  Flaw  Depth.  Failure  of  the  theory  to  predict  the 
correct  behavior  for  larger  flaws  is  attributed  to  the  snail 
flaw  approxinatlon  Inherent  in  the  nodel. 


by  the  snail  flaw  theory.  As  a  result,  it  is  believed  that  the  theory 
used  here  tends  to  overestlnate  the  effects  of  near-surface  current  per¬ 
turbations  while  underest laa ting  the  effect  of  perturbation  at  the  slot 
tip.  If  this  is  the  case,  then  a  no re  realistic  treatnent  of  the  inter¬ 
action  of  eddy  currents  with  extended  flaws  could  lead  to  flaw  response 
curves  that  change  shape  as  well  as  aagnltude  as  flaw  size  Is  Increased. 
Our  experience  with  calculations  like  those  discussed  here  leads  us  to 
believe  that  this  is  the  aost  likely  explanation  for  the  difference 
between  theory  and  experinent  shown  in  Figure  18.  Unfortunately,  the 
basic  theory  of  eddy  curent/flaw  interactions  has  not  progressed  to  the 
point  whore  aore  realistic  calculations  are  feasible. 

3.5  Slot  Tip  Effects  -  Theory 

Although  theory  cannot  provide  a  definitive  answer  concerning  the 
relative  intensity  of  near-surface  and  crack  tip  signals,  it  is  possible 
to  aake  soae  "worst  case"  and  "best  case"  estiaates  of  such  effects  based 
on  siaple  extensions  of  the  point  flaw  aodel. 
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One  such  result  is  shown  In  Figure  19.  For  this  calculation,  the 
integrated  point  flaw  model  described  in  the  preceding  section  was  used. 
This  is  a  "worst  case"  model  because  it  ignores  the  fact  that  current  is 
forced  to  flow  down  the  crack  face  away  from  the  surface,  which  would 
emphasize  crack  tip  effects  at  the  expense  of  near-surface  effects.  The 
integrated  point  flaw  model  is  therefore  expected  to  overemphasize  near¬ 
surface  perturbations  caused  by  a  crack. 


The  calculated  results  shown  in  Figure  19 
indeed  the  case.  The  curves  shown  here  are  the 
curve)  which  is  1.0  nsec  in  width  with  50  nsec 
the  flaw  response  function  (solid  curve),  both 
their  peak  values.  The  significant  point  here 
point  flaw  model  predicts  response  functions  of 
slot  depths  ranging  from  0.005  to  0.020  in.  In 
predicts  that  near-surface  effects  dominate  to 
depth  effects  cannot  be  measured. 
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input  waveform  (dotted 
rise  and  fall  times,  and 
normalized  to  unity  at 
ia  that  the  integrated 
exactly  the  same  shape  for 
other  words,  this  model 
such  an  extent  that  crack 


At  the  other  extreme,  one  might  assume  that  the  current  perturbation 
at  the  crack  tip  dominates  and  that  near-surface  effects  are  negligible. 

This  so-called  crack  tip  model  gives  the  results  shown  in  Figure  20  for 
the  same  input  pulse  used  In  Figure  19.  In  this  case,  a  very  distinct 
shift  in  decay  times  is  obtained  as  slot  depth  is  varied,  indicating  that 
depth  Information  can  be  obtained  from  an  analysis  of  pulse  decay 
characteristics. 

There  are  two  additional  points  to  be  nude  concerning  Figure  20. 

One  is  that  signal  rise  times  are  also  depth  dependent  and  that  pulse 
buildup  data  can  be  used  Just  as  well  as  decay  data  to  infer  crack  depth. 

The  second  point  is  that  decay  or  rise  time  characteristics  are  used  to 
infer  depth  rather  than  peak  arrival  times,  as  was  the  case  in  Figure  13. 
This  is  a  distinct  advantage  because  rise  and  fall  characteristics  are 
observed  over  times  of  the  order  of  hundreds  of  nanoseconds,  whereas  peak 
arrival  times  occur  in  tens  of  nanoseconds.  Thus,  for  the  particular 
pulse  shape  shown  in  Figure  20,  differences  in  decay  characteristics  are 
easily  recognized,  while  peak  arrival  time  differences  cannot  be  seen. 

The  main  question,  however,  remains  unanswered  at  this  point.  The 
integrated  point  flaw  model  used  in  Figure  19  indicates  that  depth  infor¬ 
mation  cannot  be  measured,  while  the  "best  case"  model,  the  crack  tip  cal¬ 
culation,  says  that  crack  depth  is  easily  inferred  from  the  rise  and  fall 
time  characteristics  of  the  reflected  signal.  Furthermore,  both  calcula¬ 
tions  are  suspect  because  both  are  based  on  simplified  extensions  of  the 
point  flaw  theory.  Clearly,  then,  this  situation  calls  for  an  experi¬ 
mental  study  of  signals  produced  by  flaws  of  different  depths  to  determine 
which  of  the  two  models  is  more  nearly  correct,  i.e.,  to  determine  if  flaw 
signals  are  dominated  by  near-surface  scattering  or  if  crack  tip  effects 
are  measurable. 

Before  presenting  the  experimental  results,  it  is  important  to  note 
that  experiments  with  open  slots  are  likely  to  involve  smaller  time  shifts 
than  experiments  with  closed  cracks.  The  reason  is  illustrated  in  Figure  21 
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Figure  19.  Flew  Reeponee  Based  on  an  Integrated  Point  Flaw 
Model.  The  input  current  (dashed  curve)  is  of  1  nsec  duration 
with  50  nsec  rise  and  fall  tines.  The  flaw  response  function 
(solid  curve)  is  dominated  by  near-surface  scattering  and  is 
the  sane  for  slot  depths  ranging  fron  5  to  20  nils. 
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10  15  FLAW  DEPTH  IN  MILS 


Figure  20.  Flaw  Response  Based  on  a  Crack  Tip  Model,  for 
the  Same  Input  Pulse  as  Used  in  Figure  18.  Significant 
differences  in  decay  characteristics  of  response  functions 
are  noted  for  crack  tips  at  different  depths. 


which  shows  possible  signal  propagation  paths  for  open  and  closed  flaws. 
Part  of  the  signal  from  an  open  crack  (slot)  will  arrive  at  an  earlier 
time  than  the  signal  from  a  closed  crack  of  the  same  depth  because,  in  the 
case  of  an  open  slot,  the  signal  can  return  along  a  path  that  is  mostly  in 
air.  In  the  analysis  of  experimentally  observed  time  shifts  for  open 
slots,  one  should  therefore  interpret  the  data  as  indicative  of  a  single 
propagation  path  through  the  metal  from  transmitter  to  flaw  with  a  return 
path  in  air  involving  no  time  delay.  Thus,  if  slot  tip  scattering  is 
significant,  time  shifts  should  approximate  the  time  required  for  a  signal 
to  propagate  through  the  metal  to  the  slot  tip.  If,  on  the  other  hand, 
near-surface  scattering  dominates,  then  measured  time  shifts  should  be 
much  smaller  than  the  transmitter-to-slot  tip  propagation  time. 

To  provide  a  basis  for  analyzing  experimental  data  in  this  way,  an 
analytic  solution  was  developed  for  eddy  currents  induced  by  dipole  source 
in  an  infinite  conductor.  This  solution  leads  to  the  expression 
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for  the  time  required  for  the  current  peak  to  propagate  a  distance  R  in  a 
medium  with  conductivity  o  (MKS  units). 

3.6  Slot  Tip  Effects  -  Experiment 

A  series  of  experiments  was  carried  out  to  determine  if  any  time 
resolved  flaw  information  is  available  with  the  PEC  breadboard  arrange¬ 
ment,  in  particular  to  find  out  if  any  information  from  eddy  current 
concentration  at  the  crack  tip  is  available  on  the  time  scales  implied  by 
the  above  equation.  Transmitter  probe  T4,  a  ferrite-shielded  probe  with 
only  20  turns  of  #38  gage  wire,  in  order  to  minimize  phase  distortion,  was 
used,  with  a  receiver  of  30  turns  of  If 38  gauge  wire.  The  titanium  test 
specimen  used  has  a  series  of  EDM  notches  of  constant  length  (1.25  mm)  and 
width  (0.13  mm)  but  varying  depth  from  1.25  mm  to  0.08  mm.  Since  the 
surface  area  of  each  flaw  is  the  same,  the  integrated  point  flaw  model 
predicts  no  variation  in  flaw  response.  A  pulse  length  of  approximately 
150  nsec  was  used. 

Figure  22  shows  a  sample  experiment  with  the  received  signal  in  the 
absence  of  a  flaw  and  over  the  1.25  mm  deep  flaw.  The  shapes  are  similar 
and,  in  fact,  the  peak  arrival  time  does  not  vary  substantially.  However, 
the  decay  of  the  signal  between  this  peak  and  the  first  zero-crossing  does 
exhibit  a  variation  as  the  probe  is  moved  from  unflawed  to  flawed  region, 
and  the  decay  varies  with  flaw  depth.  This  result  is  seen  more  clearly  in 
Figure  22(c)  in  which  the  no-flaw  signal  is  subtracted  from  the  flaw 
signal,  and  the  result  smoothed. 

The  difference  signals  were  smoothed  with  a  low  pass  filter  with  a 
cutoff  frequency  of  approximately  8  MHz.  The  filtered  difference  signals 
exhibit  peaks  that  decrease  in  phase  as  well  as  amplitude  as  the  flaw 
depth  decreases.  Figure  23  shows  the  peak  arrival  time  as  a  function  of 
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(a)  Signal  in  the  Absence  of  a  Flaw  (b)  Signal  in  the  Presence  of  a  Flaw 


(c)  Filtered  Difference  Signal 


Figure  22.  Received  Pulse  in  the  Absence  of  a  Flaw  (a)  Received 
Pulse  With  Probe  Positioned  to  Maximize  the  Flaw  Signal  (b) 


flaw  depth.  Also  shown  is  a  curve  of  the  time  for  an  eddy  current  pulse 
to  penetrate  from  the  surface  to  a  depth  equal  to  the  flaw  depth,  from 
Eq.  1.  The  measured  times  agree  with  the  curve  within  a  factor  of  two. 

Furthermore,  the  general  shapes  (convexity)  of  the  curves  are  simi¬ 
lar:  they  both  exhibit  time  dependence  that  is  proportional  to  the  square 
of  penetration  depth.  This  correct  dependence  confirms  that  the  times  mea¬ 
sured  are  indeed  indicative  of  eddy  current  penetration  dynamics. 

This  experiment  indicates  that  flaw  depth  information,  in  addition 
to  surface  information,  is  present  with  appropriate  signal  processing. 
Furthermore,  it  shows  behavior  that  agrees,  at  least  qualitatively,  with 
that  expected  from  the  Kahn  effect,  indicating  that  correct  modeling  of 
the  flaw  signal  requires  abandoning  the  point  flaw  model  assumptions. 
Finally,  the  experiment  and  theory  suggest  that  the  flaw  depth  information 
may  be  enhanced  by  reducing  the  direct  coupling  between  the  transmitter 
and  receiver  during  the  time  in  which  the  defect  tip  signal  reaches  the 
receiver. 

3.7  Conclusion 


Experimental  evidence  and  theoretical  calculations  have  confirmed 
that  defect  depth  information  for  surface-breaking  flaws  is  available  from 
PEC  receiver  signals.  This  Information  is  related  to  the  arrival  time  of 
the  signal  from  eddy  current  concentration  at  the  flaw  tip.  Because  the 
information  is  interpreted  from  time,  rather  than  amplitude  data,  the  mea¬ 
surement  should  be  relatively  insensitive  to  liftoff  variation. 

It  should  be  pointed  out  that  this  viewpoint  is  not  in  conflict  with 
some  previous  observations.  Waidelich  (11)  has  observed  empirically  that 
deeper  flaws  In  aluminum  influence  the  received  signal  at  later  times  than 
more  shallow  flaws.  Morris  (12)  observed,  also  empirically,  that  defect 
depth  correlated  well  with  the  observed  time  at  which  the  flaw  and  no-flaw 
signals  are  equal.  This  approach  may  be  similar  to  measuring  the  change 
in  decay  time  of  the  transmitted  pulse.  These  past  observations  have  not 
been  previously  explained  quantitatively  in  terms  of  eddy  current  paths, 
however. 

Consideration  of  this  time-domain  information  leads  to  the  question 
of  how  best  to  exploit  it  to  maximize  the  desired  signal  strength.  The 
approach  taken  in  this  project  has  been  to  take  advantage  of  the  time- 
domain  nature  of  the  signal  by  separating  in' time  the  depth  information 
from  the  surface  flaw  information  and  no-flaw  signal.  That  is,  if  signal 
transmission  could  be  halted  or  severely  diminished  before  the  flaw  depth 
signal  arrives,  then  the  slgnal-to-noise  ratio  could  be  improved. 

One  possible  current  pulse  shape  is  a  ramp,  as  shown  in  Figure  24(a). 
In  fact,  Hendrickson  and  Hansen  (13)  have  suggested  this  as  a  way  to 
generate  deeply  penetrating  eddy  currents,  then  to  receive  signals  from 
their  decay.  The  difficulty  of  this  approach,  compared  to  the  bipolar 
pulses  used  for  the  experiment  previously  described,  is  that  it  is  imprac¬ 
tical  to  stop  the  current  flow  in  the  excitor  coil  in  a  time  short  com¬ 
pared  to  the  arrival  time  expected  for  small  flaws  in  low  conductivity 
material  -  tens  of  nanoseconds  in  titanium. 
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Figure  24.  Ideal  Current  Waveforms  for  Achieving  Time-Domain 
Separation  of  Surface  Signals  and  Flaw  Depth  Information.  In 
case  (b),  a  sharp  change  in  current  is  not  required  even  though 
signal  transmission  stops  as  quickly  as  in  (a). 


The  reason  is  chat  the  coil,  if  made  capable  of  supplying  a  large 
signal,  must  have  substantial  inductance  (of  the  order  of  1  microhenry). 
The  e-folding  time  to  remove  the  current  in  an  inductor  through  a  resist¬ 
ance  is  L/R.  Therefore,  a  resistance  of  .01  ohm  is  required  to  passively 
remove  approximately  2/3  of  the  current  in  10  nanoseconds.  This  is  a 
factor  of  ten  less  than  typical  coil  resistances,  and  ignores  cable  and 
winding  capacitance.  Alternatively,  a  back  voltage  could  be  applied  to 
the  probe,  but  the  magnitude  of  this  voltage  is  inversely  proportional  to 
the  decay  time  desired.  In  other  words,  if  the  desired  decay  time  is  one 
tenth  of  the  rise  time,  then  the  decay  portion  applied  voltage  must  be 
10  times  the  rise  portion  voltage.  This  operation  mode  would  require  the 
use  of  fewer  turns  of  heavily  insulated  wire,  and  thereby  reduce  the  trans* 
mitted  power. 


Therefore,  the  selected  approach  for  this  project  is  to  use  the 
waveform  shown  in  Figure  24(b).  When  the  current  reaches  the  top  of  the 
ramp,  it  is  held  constant  for  a  time  long  enough  to  take  the  desired  mea¬ 
surements.  In  this  manner,  the  de-energizing  of  the  probe  can  take  place 
after  the  measurement  is  taken,  at  a  moderate  rate  consistent  with  its  own 
impedance.  The  voltage  induced  in  the  test  material,  however,  according 
to  Faraday's  Law,  is  proportional  to  the  rate  of  change  of  current,  and 
therefore  decays  in  the  time  that  it  takes  to  change  the  probe  current 
from  a  ramp  to  a  constant  value.  This  selected  approach  is  described  in 
the  following  section. 
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4.  ADVANCED  SYSTEM  DESIGN 


4 . 1  Breadboard  Development 

The  original  breadboard  system  was  modified  to  test  the  time  domain 
signal  processing  concepts  evaluated  in  the  latter  part  of  this  project* 

As  modified,  the  breadboard  system  consists  of  a  Texas  Instruments  Model 
6613  pulse  generator,  a  voltage-to-current  pulse  amplifier,  a  Hewlett 
Packard  Model  9826  computer,  a  Tektronix  Model  7854  digital  oscilloscope, 
and  a  Princeton  Applied  Research  Model  115  wideband  preamplifier*  In  addi¬ 
tion,  the  mechanical  scanning  system  was  modified  to  provide  digital  read¬ 
out  of  position  in  one  axis  at  approximately  0.0005  inch  resolution.  The 
modified  breadboard  system  is  shown  in  Figure  25.  In  this  system,  the 
only  components  not  commercially  available  are  the  probes  and  the  voltage- 
to-current  amplifier,  although  the  amplifier  was  constructed  out  of  commer¬ 
cially  available  discrete  components. 

The  pulse  generator  was  selected  because  it  can  provide  continuously 
variable  pulse  risetime,  falltime,  and  duration,  with  each  parameter  con¬ 
tinuously  variable  and  independently  controllable  from  10  nanoseconds  to 
at  least  1  millisecond.  The  output  is  a  voltage  pulse  of  up  to  5  volts 
peak  amplitude  into  50  ohms  resistance. 

The  pulse  generator  signal  is  amplified  by  the  voltage-to-current 
amplifier  that  was  designed  for  this  project  at  SwRI.  The  amplifier  is  a 
wideband  (approximately  50-MHz)  device  capable  of  delivering  at  least 
2.5  Amperes  into  a  load  impedance  varying  from  0  to  more  than  100  Ohms. 

This  is  achieved  at  the  cost  of  a  low  pulse  repetition  rate:  approximately 
300  pulses  per  second  at  the  maximum  current.  The  current  pulse  shape  is 
maintained  close  to  that  of  the  input  voltage  by  the  use  of  a  feedback 
circuit  that  directly  senses  the  current  delivered  to  the  probe.  This 
amplifier  is  described  in  greater  detail  and  its  circuit  schematic  is 
presented  in  Appendix  D. 

The  received  signal  is  amplified  by  a  PARC  Model  115  wideband  pre¬ 
amplifier.  The  50  Ohm  input  Impedance  of  the  preamplifier  effectively 
dampens  the  receiver  probe  response.  The  amplified  signal  is  input  to  the 
digital  oscilloscope,  which  also  serves  as  the  transient  recorder/digitizer 
for  the  system.  The  oscilloscope  is  directly  connected  to  the  HP  9826 
minicomputer  via  a  standard  General  Purpose  Interface  Bus  (GPIB).  The 
oscilloscope  waveform  acquisition  and  signal  pre-processing  can  be  con¬ 
trolled  by  the  oscilloscope  keyboard  or  the  signal  acquisition  program 
written  for  the  HP  9826.  Signal  processing  performed  in  the  oscilloscope 
Includes  multiple  signal  acquisition  and  averaging.  This  function  takes 
the  place  of  the  boxcar  averager  in  the  original  breadboard  configuration. 
Since  the  oscilloscope  reconstructs  a  high  frequency  waveform  by  repetitive 
random  sampling  of  multiple  waveforms,  it  does  not  allow  a  continous  scan 
to  be  made  of  the  signal  at  one  point  in  time.  Hardcopy  output  is  provided 
by  a  thermal  printer  attached  to  the  HP  9826. 
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Figure  25.  The  Modified  Pulsed  Eddy  Current  Test  System.  This 
breadboard  system  has  been  constructed  to  allow  use  of  monopolar 
pulses  with  variable  rise/fall  times  and  duration,  and  more  flex¬ 
ible  (software  controlled)  feature  extraction.  Only  the  SwRI- 
constructed  voltage-to-current  pulse  amplifier  is  not  commercially 
available. 


This  breadboard  system  provides  greater  flexibility  in  pulse  shaping 
and  signal  processing  at  the  cost  of  reduced  probe  current,  reduced  pulse 
repetition  rates,  and  inability  to  provide  real  time  scans  of  flaws  speci¬ 
mens.  The  limited  peak  current  of  2.3  amperes  compares  to  more  than 
7  amperes  typically  available  with  probes  connected  to  the  Velonex  pulser, 
where  the  current  was  effectively  limited  only  by  the  need  to  maintain  the 
peak  voltage  below  the  arcing  point  of  the  probe  windings.  Thus,  the  flaw 
detection  capability  was  sacrificed  and  cannot  be  directly  compared  with 
that  of  the  original  pulsed  eddy  current  breadboard  system  or  the  conven¬ 
tional  eddy  current  system. 

Figure  26  shows  the  pulse  quality  available  through  a  one-half  turn 
coil  of  wire,  essentially  a  short  circuit.  The  bottom  trace  is  the  input 
voltage  waveform,  and  the  top  trace  is  the  output  current,  at  one  ampere 
per  division.  The  horizontal  scale  is  5U0  nanoseconds  per  division. 

4 . 2  Probe  Development 

Transmitter  probe  response  to  the  pulsing  system  was  studied.  There 
were  two  areas  of  concern.  One  is  that  the  voltage  swings  necessary  to 
drive  the  probe  at  the  desired  waveform  be  kept  below  the  limited  capabil¬ 
ity  of  the  pulser  (less  than  200  volts).  Second  is  the  desire  to  make  the 
probe  coil  current  follow  the  pulser  output  current.  These  two  currents 
may  not  be  equal  because  capacitance  in  the  probe  cable  and  coil  Itself 
will  shunt  part  of  the  current.  It  has  been  experimentally  observed  that 
a  lower  resonant  frequency  causes  greater  output  voltage  swings  and  thus 
reduces  the  ability  of  the  feedback  circuit  to  track  the  input  current. 
Figures  27(a)  and  27(b)  show  this  effect  for  a  b-MHz  resonant  frequency 
probe  (Figure  27(a))  and  a  14-MHz  probe  [Figure  27(b)). 


Figure  26.  Input  Voltage  Waveform  (Bottom  Trace)  and  Output  Current 
Waveform  (1  Ampere  per  Division)  of  the  Voltage-to-Current  Amplifier. 
Output  is  through  a  one-half  turn,  one-half  inch  long  coil  of  wire. 
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Figure  27.  Probe  Current  (Upper  Trace)  Follows  Input  Voltage 
(Lower  Trace)  Better  When  a  Probe  With  a  Higher  Resonant 
Frequency  [(b),  14.3  MHz]  Probe  is  Used  than  When  a  Lower 
Resonant  Frequency  Probe  ((a),  6  MHz]  is  used.  Probes  have 
sistflar  construction,  but  use  different  cables.  The  pulse 
rise  tiae  is  400  nanoseconds. 
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In  order  to  study  the  coil  response  to  an  ideal  probe  input  wave- 
form,  the  probe  was  modeled  as  a  capacitance  in  series  with  an  inductance. 
The  result  for  the  difference  between  current  into  the  probe  and  through 
the  coil  is 


I  -  I0  -  sin  apt 
usqT 


sin  apT  cos  u>0(t-T) 
|*>0T 


+  1  -  cos  mqT  S£n  Wq  (t-T) 
“0T 


where  u>q  is  the  resonant  frequency  of  the  probe  in  radians  per  second  [u>q2 
-  1/(LC) ] ,  T  is  the  duration  of  the  ramp,  t  is  time  since  beginning  of  the 

pulse,  and  IQ  is  the  current  in  amperes  at  the  end  of  the  ramp  portion  of 

the  pulse  and  during  the  constant  portion. 

This  analysis  Indicates  that  in  order  for  the  coll  current  to  accu¬ 
rately  track  the  pulser  output,  the  following  condition  must  be  met: 

fT  »>i 
2 

where  f  “  Wq/2*  is  the  resonant  frequency  in  Hertz. 

An  effective  way  to  increase  the  probe  resonant  frequency  (and  band¬ 
width)  without  decreasing  its  transmitting  power  is  to  reduce  the  capaci¬ 
tance  of  the  pulser-to-coil  cable.  Figure  28  shows  the  effect  of  cable 
capacitance  for  probe  T4  (20  turns  of  #38  gage  wire  on  a  ferrite  core,  fer¬ 
rite  shielded).  Without  a  cable,  the  apparent  Inductance  of  the  probe  is 

almost  constant  to  15  MHz.  With  standard  50-ohm  impedance  RG/U  174  cable, 
the  resonant  frequency  is  lowered  to  10  MHz.  This  cable  has  a  nominal 
capacitance  of  30.8  pf  per  foot.  With  RG/U  62B  cable  at  13  pf  per  foot, 
the  resonant  frequency  is  Increased  to  14  MHz.  Finally,  with  the  use  of  a 
low  capacitance  cable  (RG/U-114A)  rated  at  10  pf  per  foot,  the  probe  has  a 
resonance  at  17  MHz.  It  was  decided  to  use  RG/U  62B  cable  for  probe  devel¬ 
opment  efforts  because  the  marginal  Improvement  of  the  R6/U-114  cable  is 
offset  by  greater  bulk,  stiffness,  and  electrical  noise. 

Figure  29  shows  the  current  waveform  that  is  achievable  with  the 
10-Mllz  probe,  using  a  250-nanosecond  ramp,  and  the  resulting  power  spec¬ 
trum.  It  is  seen  that  the  input  power  Intensity  is  negligible  for  radial 
frequencies  above  n/T.  Figure  30  shows  in  comparison  the  ideal  frequency 
magnitude  spectrum.  The  magnitude  peaks  above  the  first  zero  actually 
represent  only  1  percent  of  peak  power.  Therefore,  it  is  concluded  that 
probes  can  be  constructed  to  allow  the  pulsing  unit  to  deliver  a  reasonably 
accurate  waveform  to  the  transmitter  probe  coil. 
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Figure  28.  The  Variation  of  Observed  Inductance  of  a  Single  Transmitter  Coil  as  a 
Function  of  Frequency,  When  Measured  With  Short  Leads  and  With  2  Feet  of  Different 
Coaxial  Cables.  The  cable  substantially  decreases  probe  bandwidth,  but  this  effect 
may  be  minimized  by  use  of  low-capacitance  cable  (RG-62B  or  RG-114A). 
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Figure  30.  Magnitude  of  the  Fourier  Transform  of  the  Ideal 
Transmitter  Current  Waveform  T  is  the  Duration  of  the  Ramp 
Portion.  Note  that  at  frequencies  above  the  first  zero,  the 
peak  magnitude  is  10%  of  the  first  peak,  so  the  energy  is 
only  1  percent. 
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4.3  Evaluation  Results 


Tests  were  conducted  with  the  new  pulsing  system  and  transmitter 
probe  T4  modified  with  the  low  capacitance  coaxial  cable  and  the  addition 
of  an  absolute  receiver  coll  placed  parallel  to  the  transmitter.  As  in 
previous  experiments,  the  receiver  coll  was  aligned  first  visually  and 
then  by  minimizing  the  transmitter-to-receiver  coupling  when  placed  over 
an  unflawed  surface. 

A  sample  of  test  results  is  shown  in  Figures  31  and  32.  Figure  31 
shows  the  receiver  signal  in  the  absence  of  a  flaw,  during  the  ramp  and 
constant  value  portions  of  the  transmitted  signal.  Figure  32  waveforms 
are  the  difference  between  a  flaw  signal  and  the  no-flaw  signal  of  Fig¬ 
ure  31.  This  approach  is  similar  to  that  of  Section  3;  however  in  this 
case  the  resulting  curve  is  not  smoothed.  It  is  seen  that  with  the  use  of 
this  new  pulse  shape,  the  difference  in  decay  characteristics  predicted  by 
the  crack  tip  model  as  shown  in  Figure  20  can  now  be  directly  observed. 
Figure  33  shows  the  inverse  of  the  slope  of  the  decay  portion  of  the  curve 
as  a  function  of  flaw  depth.  The  inverse  slope  is  used  because  this  quan¬ 
tity  is  a  measure  of  the  time  delay  of  the  signal  in  reaching  any  given 
value,  and  should  therefore  increase  with  flaw  depth.  It  is  seen  that 
flaw  depth  corresponds  to  this  measured  parameter,  in  agreement  with  the 
crack-tip  model  calculations  shown  in  Figure  20.  A  more  quantitative 
comparison  is  difficult  to  make  because  of  the  variation  between  ideal  and 
achieved  pulse  shapes. 

Two  encouraging  observations  can  be  made  from  these  results.  First, 
the  basic  predictions  made  from  the  crack  tip  model  have  been  confirmed. 
Second,  this  result  has  suggested  a  way  to  measure  flaw  depth  that  appears 
to  be  more  sensitive  than  that  used  in  Section  3.  In  this  approach,  the 
smallest  target  flaw  appears  on  the  scale,  with  a  finite  time,  in  contrast 
to  the  result  obtained  with  the  initial  breadboard  system,  in  which  the 
measured  time  delay  of  the  0.005-inch  deep  flaw  was  0. 


Figure  31.  The  No-Flaw  Receiver  Signal  Used  in  a  Series  of 
Flaw  Depth  Measurement  Tests.  The  amplitude  drops  sharply 
when  the  transmitted  current  becomes  constant. 


Figure  32.  Receiver  Flaw  Minus  No-Flaw  Signal,  for  Different  Depth 
Notches.  The  slope  of  the  decaying  portion  of  the  signal  decreases 
with  increasing  flaw  depth,  as  predicted  by  the  crack-tip  model 
results  of  Figure  20. 
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Figure  33.  Inverse  of  the  Slope  of  the  Decay  Portion  of  Flaw 
waveform  shown  In  Figure  32.  The  curve  Is  siailar  to  that 
obtained  from  the  bipolar  pulse  shown  In  Figure  23  confining 
the  presence  of  depth  information  on  an  eddy  current  diffusion 
tine  scale. 


5.  COMPARISON  WITH  C-W  RESULTS 


A  conventional  eddy  current  analysis  was  carried  out  on  the  flaw  samples 
in  order  to  provide  a  comparison  between  current  state-of-the-art  practice 
and  PEC  testing  capability  as  developed  in  this  project.  The  eddy  current 
instrument  used  was  a  Nortec  NDT-19.  This  instrument  is  a  single  fre¬ 
quency  device,  with  a  maximum  frequency  of  A. 6  MHz.  The  probe  used  was  a 
Nortec  SP-5000.  This  probe  is  a  conventional  absolute  probe  with  windings 
parallel  to  the  surface.  Throughout  the  experiments,  the  auto  balance 
mode  of  the  NDT-19  was  used  to  null  the  probe. 

The  test  fixture  of  the  modified  PEC  breadboard  system  was  used  to  hold 
the  probe.  During  initial  experiments,  it  was  found  that  surface  varia¬ 
tion  in  the  test  specimens  created  substantial  liftoff  signal  variation, 
so  the  probe  fixture  was  modified  to  maintain  constant  probe  contact  with 
the  surface  by  spring-loading  the  probe. 

Tests  were  conducted  by  performing  a  series  of  scans  in  the  Y-direction, 
that  is,  parallel  to  the  length  of  the  specimen  blocks.  At  the  end  of  each 
Y-scan,  the  probe  position  was  incremented  0.025-inch  in  the  X-direction 
before  another  Y-scan  was  conducted.  Limit  switches  on  the  test  fixture 
provided  signals  to  the  data  acquisition  equipment  to  ensure  that  each 
Y-scan  started  and  ended  at  the  same  value  of  X. 

The  eddy  current  instrument  was  adjusted  such  that  the  initial  liftoff 
signal  was  in  the  horizontal  left,  or  -H,  direction.  This  allowed  the 
vertical  component  of  the  signal  to  be  used  as  the  indicator  of  flaw 
characteristics.  This  approach  makes  the  vertical  signal  more  sensitive 
to  flaws.  Both  horizontal  and  vertical  components  of  the  instrument 
Impedance  plane  display  were  sent  to  the  data  acquisition  system.  Data 
acquisition  and  processing  were  conducted  by  the  use  of  a  Biomat ion 
digitizer  connected  to  a  Data  General  Nova  series  computer. 

Several  different  presentation  formats  were  developed.  Figure  34  shows 
two  representations  of  flaws  A  through  D  on  the  Ti-6-4  surface  flaw  detec¬ 
tion  standard.  This  figure  was  produced  by  a  single  scan  in  the  Y-direction 
through  the  center  of  each  flaw.  The  flaws  are  all  0.010-inch  long  by 
0.005-inch  wide.  The  depths  of  flaw  A  through  D  are,  respectively,  0.050, 
0.035,  0.020,  and  0.005-inch.  Flaw  E,  which  is  0.003-inch  deep,  was  not 
detectable.  Figure  34(a)  is  a  re-creation  of  the  impedance  plane  display 
of  the  eddy  current  instrument...  Figure  34(b)  shows  the  separate  horizontal 
and  vertical  components  of  the  impedance  plane  display  as  a  function  of 
scan  position.  In  both  presentation  forms,  it  is  seen  that  the  response 
of  flaws  A,  B,  and  C  is  very  similar. 

Figure  35  shows  a  plan  view  created  by  presenting  the  vertical  component 
of  the  impedance  plane  for  a  series  of  scans  over  the  same  flaws  as  shown 
in  Figure  34.  It  should  be  noted  that  each  plan  view  in  the  figure  was 
automatically  scaled  by  a  different  factor  in  order  to  maximize  resolution 
in  the  dlgital-to-analog  converter.  The  actual  signal  level  represented 
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Figure  34.  Conventional  Eddy  Current  Scans  of  Flaws  A  Through  D  on 
the  Surface  Flaw  Detection  Standard  (E  is  not  Detectable).  Top  figure 
is  computer  reconstruction  of  eddy  current  instrument  Impedance  plane 
display.  Bottom  figure  is  separate  plots  of  vertical  and  horizontal 
impedance  plane  components.  The  axes  are  adjusted  such  that  liftoff 
is  in  the  negative  horizontal  direction.  The  flaw  signature  is  similar 
for  flaws  A,  B,  and  C,  which  have  equal  surface  areas,  but  depths 
varying  from  0.050  to  0.020  inch. 
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Figure  35.  Plan  View  Representation  of  the  Vertical  (Mimimum  Liftoff 
Effect)  Component  of  the  Impedance  Plane  Display  for  the  Same  Flaws  as 
In  Figure  34.  Scans  are  spaced  .025-inch  apart.  Each  scan  set  is 
scaled  by  the  maximum  signal  in  the  set,  as  shown  at  the  lower  left  of 
each  view. 


in  each  view  is  proportional  to  the  maximum  signal  level  in  that  view, 
shown  at  the  bottom  of  each  scan  view.  From  these  scans,  it  is  seen  that 
the  maximum  flaw  signal  occurs  when  the  probe  is  centered  over  the  end  of 
the  flaw.  When  the  probe  is  centered  over  the  middle  of  the  flaw,  the 
eddy  currents  are  not  as  disturbed,  because  the  probe  is  larger  than  the 
flaws,  and  the  eddy  current  density  induced  below  the  center  of  the  probe 
is  small. 

Figure  36  shows  the  same  plan  view  display  as  described  above,  for  a 
series  of  flaws  (G,  H,  I,  and  J)  with  a  constant  aspect  ratio  of  2  to  1 
(length  to  depth).  Flaw  F  in  this  series  is  the  same  size  as  Flaw  E  in 
the  previous  series  and  is  also  not  consistently  detectable.  It  is  seen 
that  the  variation  in  response  is  much  more  dramatic  than  for  the  previous 
series  ABCD,  which  consists  of  flaws  with  equal  surface  area.  This  point 
is  illustrated  in  Figure  37,  which  shows  the  maximum  response  as  a  func¬ 
tion  of  flaw  depth  for  each  series. 

Figure  38  shows  data  from  the  probe  system  being  developed  for  the  RFC 
program,  for  the  same  set  of  flaws  as  in  Figure  36.  A  conventional  2  MHz 
eddy  current  coil  is  used,  but  it  is  mounted  in  an  air-bearing  probe  to 
provide  a  relatively  constant  liftoff.  The  results  are  seen  to  be  con¬ 
sistent  with  those  described  above.  In  particular.  Flaws  G,  H,  I,  and  J 
are  easily  detectable  with  increasing  signal  amplitudes,  but  Flaw  F  is  not 
clearly  distinguishable  from  noise.  * 

Two  significant  observations  can  be  made  from  this  conventional  data. 
First,  the  conventional  system  set  up  in  this  project  is  capable  of 
meeting  the  target-size  flaw  detection  goals  set  up  for  the  PEC  system, 
and  shows  flaw  detection  capability  similar  to  that  of  the  RFC  program 
eddy  current  system. 

Second,  the  response  of  the  conventional  system  is  approximately  the  same 
for  flaws  of  equal  surface  area,  for  depths  ranging  from  0.020-lnch 
(0.3  mm)  to  0.050-inch  (1.25  mm),  and  varies  substantially  with  flaw 
length.  Therefore,  it  is  expected  that  conventional  continuous-wave  eddy 
current  inspection  could  be  developed  to  provide  an  indication  of  surface- 
breaking  flaw  surface  area,  but  it  does  not  appear  as  promising  for  flaw 
depth  characterization. 


Figure  36.  Plan  View  Representation  of  the  Vertical  (Minimus  Liftoff 
Effect)  Component  of  the  Impedance  Plane  Display  for  Flaws  G,  H,  1, 
and  J.  Each  scan  set  is  scaled  by  the  maximum  signal  in  the  set  as 
shown  at  the  lower  left  of  each  view.  There  is  a  much  greater  varia¬ 
tion  in  signal  response  than  for  the  equal  surface  area  flaws  A,  B,  C 


6.  CONCLUSIONS  AND  RECOMMENDATIONS 


It  has  been  shown  that  pulsed  eddy  current  examination  is  capable  of 
detecting  the  target  flaws  in  low  conductivity  alloys,  and  that  it  offers 
a  promising  way  to  estimate  the  depth  of  surface-breaking  flaws* 

It  has  been  successfully  demonstrated  that  the  pulsed  eddy  current  bread¬ 
board  system  is  capable  of  detecting  flaws  as  small  as  0.010  inch  by 
0.005  inch  in  a  titanium-6-4  surface.  The  equipment  required  to  set  up 
such  a  PEC  system  is  commercially  available. 

By  using  the  unique  time-domain  characteristics  of  PEC  signals,  the  corre¬ 
lation  receiver  flaw  signal  decay  with  flaw  depth  was  shown.  This  is  pos¬ 
sible  because  the  eddy  current  density  is  concentrated  at  the  tip  and 
depleted  at  the  surface  of  a  surface-breaking  flaw.  Theoretical  calcula¬ 
tions,  when  revised  to  Include  this  Kahn  effect,  showed  reasonable  agree¬ 
ment  with  the  measured  data  for  flaw  depth. 

Because  the  concept  of  the  depth  measurement  technique  is  relatively  unex¬ 
plored,  it  is  recommended  that  additional  research  be  conducted  in  this 
area.  Several  areas  of  improvement  need  to  be  investigated.  First,  the 
breadboard  pulsing  unit  containing  the  voltage-to-current  amplifier  needs 
to  deliver  larger  currents  than  is  now  possible  at  equal  or  better  band¬ 
width.  This  would  make  possible  a  comparison  between  the  flaw  detection 
capability  of  the  ramp  pulse  system  and  the  more  conventional  bipolar 
pulse  system  that  was  set  up  at  the  beginning  of  the  project. 

Second,  theory  development  should  be  pursued,  in  the  area  of  flaw  response. 
The  models  that  were  developed  did  successfully  predict  the  spatial  and 
temporal  variation  of  transmitted  magnetic  fields  due  to  shielded  eddy 
current  probes,  but  the  integrated  point  flaw  model  is  clearly  inadequate 
to  describe  flaw  signal  behavior  when  flaw  characterization  is  the  goal. 
Even  though  the  eddy  current  paths  of  interest  are  three-dimensional  in 
nature,  the  authors  believe  that  an  appropriate  step  is  to  develop  an 
analytical  two-dimensional  model  that  could  simulate  a  long  crack  with  a 
finite  depth.  This  is  the  simplest  model  that  would  self-consistently 
Include  the  Kahn  effect.  Theory  development  in  this  area  may  help  provide 
guidance  to  further  probe  development  for  depth  measurement,  including  the 
use  of  differential  receiver  coils. 

Third,  the  depth  measurement  approach  needs  to  be  tried  out  on  naturally- 
occurring  flaws.  This  test  is  Important,  because  it  is  theorized  that 
this  measurement  depends  on  the  flaw  geometry,  and  not  only  the  flaw  sur¬ 
face  area,  as  for  example  appears  to  be  the  case  with  electric  current 
perturbation  and,  to  a  lesser  extent,  conventional  eddy  current 
examination. 
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APPENDIX  A 


ANALYTICAL  MODELS 


Development  of  the  Simplified  Model 


The  purpose  of  developing  the  simplified  analytical  model  was  to  obtain 
approximate  mathematical  solutions  that  describe  the  effect  of  the  shield 
on  the  transient  field  produced  by  an  arbitrary  excitation  pulse*  Of  pri¬ 
mary  Interest  is  the  focusing  effect,  i.e.,  the  narrowing  of  the  spatial 
distribution  of  field  and  the  shielding  of  the  detector  from  the  excitation 
pulse.  To  achieve  the  stated  purpose,  the  field  reflected  by  the  specimen 
was  ignored,  and  the  model  employed  was  a  cylindrical  excitation  coll  of 
finite  length  surrounded  by  a  cylindrical  shield  of  the  same  length.  The 
probe  configuration  of  the  excitation  coil  used  Is  as  shown  in  Figure  39. 

An  approximate  analytical  solution  was  obtained  by  first  solving  for  the 
transient  current  induced  in  a  shield  of  infinite  length,  and  then  substi¬ 
tuting  this  current  in  the  integral  expression  for  the  vector  potential 
for  a  shield  and  coil  of  finite  length.  The  solution  thus  obtained  was 
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where  Bz  is  the  vertical  component  of  the  flux  density  (in  the  sense  of 
Figure  39),  JQ  is  the  excitation  current,  is  the  coil  radius,  P2  is  the 
inner  radius  of  the  shield,  and  G  is  a  rather  complicated  geometrical  fac¬ 
tor  involving  elliptic  integrals.  (8)  This  expression  applies  to  thin 
shields  such  that  T  <<  |x|,  where  T  is  the  shield  thickness.  It  also 
applies  only  to  shielding  of  the  excitation  pulse;  the  effect  of  the  field 
that  diffuses  through  the  shield  has  been  ignored.  Use  of  this  result  is, 
therefore,  limited  to  studies  of  the  spatial  distribution  of  the  primary 
field. 

To  estimate  the  time-dependent  field  that  diffuses  through  the  shield  to 
the  detector  position,  the  solution  for  an  infinitely  long  shield  was 
employed  with  the  additional  assumption  that  the  effect  of  shield  curva¬ 
ture  could  be  ignored.  This  led  to  an  exponentially  convergent  series 
solution  for  the  field  at  the  outside  surface  of  the  shield. 

These  approximate  analytical  solutions  were  used  to  study  the  effects  of 
probe  geometry  on  the  focusing  of  the  primary  field  and  the  time-dependence 
of  the  field  at  the  outside  surface  of  the  shield.  A  typical  result  illus¬ 
trating  the  focusing  effect  is  shown  in  Figure  40.  This  is  a  plot  of  the 
geometrical  factor  P2G(xiPi)  -  p^G(x^P2)  as  a  function  of  distance  from 
the  center  of  the  coll  along  a  line  0.025  mm  from  the  end  of  the  probe. 

The  sharp  negative  peak  at  about  1.8  mm  occurs  in  the  region  between  the 
coil  and  shield. 


SHIELD 


Figure  39.  Simplified  Model  of  Pulsed  Eddy  Current  Probe. 
This  model  was  used  to  evaluate  the  field  in  front  of  and 
alongside  the  probe.  The  finite  current  coil  is  modeled  as 
a  thin  current  sheet  at  pi. 
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The  time-dependent  field  just  inside  and  outside  a  shield  of  0.25-mm  thick¬ 
ness  is  shown  in  Figure  41.  For  a  shield  this  thin,  the  calculation  shows 
that  the  diffusion  component  of  the  field  just  outside  a  shield  peaks  about 
1  usee  after  the  initial  pulse.  Other  calculations  indicate  that  increas¬ 
ing  the  shield  thickness  to  1.5  mm  causes  the  diffusion  time  to  increase 
to  about  9  usee  and  reduces  the  amplitude  of  the  diffusion  pulse  by  about 
a  factor  of  200. 

These  results  indicate  that  copper  shielding  provides  both  focusing  of  the 
excitation  field  and  shielding  of  receiver  coils  from  coupling  directly 
with  the  excitation  field. 

Development  of  Detailed  Model 


The  detailed  model  is  a  small  flaw  theory  in  the  sense  that  the  current 
perturbation  caused  by  the  flaw  is  approximated  by  a  current  dipole  in 
exactly  the  same  way  as  in  the  theory  of  Burrows  (9J  and  Dodd  et  al.  (10) . 

The  major  difference  between  this  new  small  flaw  model  and  the  Burrows' 
theory  is  that  we  start  with  the  Fourier  or  Laplace  transform  of  the 
general  time-dependent  form  of  Maxwell's  equations  rather  than  the  special 
case  of  a  sinusoidal  driving  current. 

Application  of  the  reciprocity  theorem  and  the  introduction  of  the  small 
flaw  scattering  matrix  5  then  follow  as  in  Burrows'  dissertation,  leading 
to  the  flaw  Impulse  response  function  given  by 
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where  a  is  the  conductivity,  Vf  is  the  flaw  volume,  a  is  the  flaw  scatter¬ 
ing  matrix,  and  the  E/1  functions  are  electric  fields  at  the  flaw  position 
caused  by  current  pulses  of  Infinitesimal  duration  in  the  transmitter  (t) 
and  receiver  (r).  For  an  excitation  current  I(t)  in  the  transmitter,  the 
voltage  pulse  in  the  receiver  is  given  by  the  convolution  of  I(t)  with 
Av(t),  i.e., 
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The  first  step  in  a  prediction  of  flaw  response  is  the  calculation  of  the 
vector  potentials  A0  and  radial  flux  densities  3A0/3z  in  air  for  an  axisym- 
metric  transmitter  and  an  axisymmetric  differential  receiver.  Input  vari¬ 
ables  for  these  calculations  are  coil  length,  radius  of  the  permeable  core. 
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Figure  41.  Time-Dependent  Field  Alongside  a  Probe  With  a  0.25-mm  Thick. 
Shield.  It  takes  approximately  one  microsecond  after  the  initial  pulse 
for  the  field  peak  to  diffuse  through  the  shield.  The  current  sheet  is 
1.6  mm,  and  the  shield  is  1.7  mm  from  the  probe  centerline. 
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winding  thickness  and  shield  radius.  For  a  differential  coil,  the  winding 
consists  of  any  number  of  alternating  layers  of  positive  and  negative  cur¬ 
rent.  The  theory  used  for  these  calculations  was  the  same  as  that  used  in 
the  simplified  model. 

Once  the  fields  in  air  are  known,  the  next  step  is  the  calculation  of  E/1 
functions  for  the  transmitter  and  receiver.  This  may  be  done  by  direct 
numerical  integration  of  the  following  exact  expression: 
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where  AQ  and  its  derivative  are  evaluated  at  P  '  and  the  Fn  are  Integrals 
of  ordinary  Bessel  functions.  Two  of  these  Integrals,  and  F3,  can  be 
evaluated  analytically,  but  the  third  requires  numerical  integration.  The 
exact  calculation  of  E/I  therefore  requires  rather  extensive  numerical 
computation. 

In  numerical  studies  of  the  various  quantities  Involved  in  this  integral, 
it  was  discovered  that  the  Fn  functions  are  strongly  peaked  at  p''  »  p. 

This  finding  then  led  to  the  following  approximate  result  in  which  AQ  and 
9A0/3z  are  evaluated  at  p: 
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where  IQ  and  1^  are  modified  Bessel  functions.  Numerical  tests,  one  of 
which  is  shown  in  Figure  42,  show  that  the  approximate  formula  gives 
results  that  are  within  about  20  percent  of  the  results  obtained  by  numer¬ 
ical  integration  except  for  small  p  where  E/I  itself  is  small.  Because 
use  of  the  approximate  formula  reduces  computer  time  by  about  a  factor  of 
1000,  all  parametric  flaw  response  studies  were  based  on  the  approximate 
calculation  of  E/I. 

The  next  step  is  the  calculation  of  Av(t)  according  to  equation  (A-l). 

The  time  integral  is  evaluated  numerically  and  the  results  are  retained  as 
a  disc  file.  The  final  step  is  then  given  by  equation  (A-2),  which,  again, 
is  evaluated  numerically. 


Figure  42.  Comparison  of  the  Analytic  Approximation  to  the 
Exact  Expression  of  Current  Density 
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APPENDIX  B 

AN  ALTERNATIVE  APPROACH  FOR  COMPUTING  EDDY  CURRENTS: 
CASE  OF  THE  DOUBLE-LAYERED  PLATE 


An  Alternative  Numerical  Approach  for  Computing 
Eddy  Currents:  Case  of  the  Double-Layered  Plate 

MARTIN  J.  SABLIK,  member,  ieee.  ROBERT  E.  BEISSNER,  and  ATHENE  CHOY 


Abstract- A  numerical  approach  ia  presented  for  caaaa  where  induced 
eddy  current  denritfec  are  to  be  computed  in  axiiymmetric  situations. 
Conductors  are  represented  as  stacks  of  inductively  coupled,  concen¬ 
tric,  metal  rings  in  which  current  to  excited  by  an  external  cod.  Matrix 
techniques  ate  used  to  solve  the  resulting  aet  of  circuit  equations.  Sinu¬ 
soidally  induced  currents  in  a  sandwich  plate  consisting  of  two  layers  of 
different  metals  are  discussed  as  an  application.  Comparison  to  made 
with  a  finite  element  calculation  and  with  results  obtained  via  an  ana¬ 
lytical  approach.  Non  sinusoidal  excitation  to  also  briefly  discussed. 
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I.  Introduction 

N  THIS  PAPER  we  discuss  a  numerical  approach  to  axisym- 
metric  problems  involving  the  computation  of  eddy  current 
distributions  induced  in  conductors  by  a  driving  current  in  c 
cylindrical  coil.  Each  conductor  is  represented  by  stacks  of 
inductively  coupled  metal  rings  concentric  with  the  driving 
coil.  Matrix  techniques  are  employed  in  solving  the  set  of  cir¬ 
cuit  equations  used  to  represent  the  interaction  between  rings 
in  the  conductors  and  the  external  driving  current.  This  “cou¬ 
pled  circuit  approach”  is  a  generalization  of  earlier  work  by 
Chang  and  Smith  (1] ,  (2) . 

The  paper  is  organized  as  follows.  Section  II  discusses  the 
coupled  circuit  approach.  Section  HI  presents  numerical  re¬ 
sults  for  the  case  of  eddy  currents  induced  in  a  sandwich  plate 
consisting  of  two  layers  of  different  metals  and  compares  these 
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Fig.  1.  Example  of  axisymmctric  configuration  which  can  be  treated 
by  coupled  circuit  approach.  Metal  plate  and  shield  arc  subdivided 
into  concentric  sucks  of  rings,  as  shown  in  this  two-dimensional 
cron  section. 


results  to  those  obtained  by  Palantsamy  et  at.  (3] ,  (4]  via  the 
finite  element  method  [5) ,  [6] .  Section  IV  compares  the  ana¬ 
lytic  solution  of  Dodd  eial.  [7]  -|9]  to  the  coupled  circuit  re¬ 
sults  for  the  sandwich  plate.  Section  V  briefly  discusses  exten¬ 
sion  of  the  coupled  circuit  approach  to  nonsinusoidal  driving 
currents. 


II.  The  Coupleu  Circuit  Approach 
When  current  flows  in  a  circular  coil,  the  vector  potential 
produced  by  it  has  only  one  component  A*  in  cylindrical  co¬ 
ordinates.  Thus  currents  induced  in  a  conductor  tend  to  flow 
in  circular  paths  concentric  with  the  axis  of  the  excitation  coil. 
The  same  situation  prevails  when  the  excitation  coil  is  a  multi¬ 
turn  tightly  wound  cylindrical  coil.  Thus,  when  conductors  are 
positioned  axisymmetrically  with  respect  to  the  excitation 
coil,  the  conductors  can  be  modeled  as  stacks  of  inductively 
coupled,  concentric,  conducting  rings. 

fig.  1  shows  a  typical  problem.  A  coil  is  surrounded  by  a 
conducting  shield  and  positioned  over  a  conducting  plate.  The 
plate  is  modeled  as  a  disk  broken  horizontally  and  vertically 
into  stacks  of  concentric  conducting  rings.  In  Chang  and 
Smith's  model  (I|  of  a  conducting  plate,  a  uny  cylindrical 
hole  was  inserted  in  the  center  of  the  plate.  This  hole  is  not 
necessary  because  the  vector  potential jjoes  to  zero  along  the 
central  axis  (lOj  and  hence  so  does  E  -  -dA/dt  and  ./  *  at. 
Our  computer  calculations  produce  essentially  equivalent  re¬ 
sults  with  and  without  the  hole. 

The  problem  depicted  in  Fig.  I  is  of  greater  complexity  than 
that  treated  by  Chang  and  Smith  [I],  [2]  in  that  I)  a  multi¬ 
turn  excitation  coil  is  considered,  2)  i  vertical  distribution  of 
eddy  currents  in  the  plate  is  obtainable,  and  3)  the  effect  of  a 
conducting  shield  about  the  driving  coil  can  be  considered. 


Eddy  current  distributions  in  both  the  plate  and  the  shield  are 
obtainable.  In  another  type  of  general  axisymmetric  problem, 
the  shield  can  be  replaced  by  a  secondary  coil  and  currents  in 
the  secondary  can  be  obtained.  The  plate  can  also  be  divided 
into  layers  of  different  conducting  media. 

Each  of  the  rings  is  replaced  by  its  corresponding  circuit  rep¬ 
resentation.  The  ith  ring  is  assigned  a  self-inductance  L,  and  a 
resistance  Rf.  Mutual  inductance  Ml0  exists  between  the  r'th 
ring  and  the  excitation  coil.  In  addition  the  tth  ring  is  coupled 
to  all  other  rings  via  mutual  inductances  My  =  M/t,  where  i* 
1,2  •  ■  ■  N  and  i  *  j  and  where  there  are  N  rings  in  total. 

The  currents  in  the  rings  are  then  completely  determined 
[  1  ]  by  the  following  sets  of  equations  in  phasor  form  for  sinu¬ 
soidal  angular  frequency  co: 

(R  |/|  +/wf.|/,)  +  /wM12/2  +•  •  ■+/cuM,w/yv 

=  */uM  |q/q 

/cuAf2I/|  +  (f?2/2  +/a>Z.2/2)  +  •  ■  •  +/o oM1NlN 

*  -/w/tfjo/o 

/uiMfj  |  / 1  +/uAf/v2/2  ■*•  ••  +  {RnIn  +Ao£/v//v) 

*  (1) 
where  /0  is  the  current  in  the  excitation  coil.  After  dividing 
both  sides  of  these  equations  by  /to,  the  equations  may  be 
written  in  matrix  form  as 


mi/i— 


where 

m 


L\  -  iRil<t)  Af|  2 

Mu  Li  -  /Ri/uj 

Mn  i  Ms 2 


Mis 

Mis 


Ls  '  iRnIv I 


(2) 


(3) 


and  [/]  is  the  column  matrix  [/|,  /2,  •  ■  •  ,/s]  and  [V]  the 
column  matrix  [Afto/o,  Af20/o.  ■  ■  •  ,MN0f, 0) .  Equation  (2) 
can  then  be  solved  as 


[/]■-  win  (4) 

To  calculate  the  mutual  inductance  between  two  circular 
filament  conductors  with  the  same  axis,  the  following  result  is 
used  [11),  (12): 

Mi  2  *  Mo  [(I  -  *)  K(k)  -  -  £(*)]  ,  (5) 


where  r,  and  r2  are  the  radii  of  the  circular  filaments,  K(k) 
and  E(k)  are  complete  elliptic  integrals  of  the  first  and  second 
kind,  and 


where 


P1  ~  (ri  -  M2  +  sJ.  (7a) 

<12  =(rj  +r>)2  +  JJ.  (7b) 

and  s  is  the  distance  between  the  planes  of  the  filaments. 
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Following  Rota  and  Grover  [12] ,  Chang  [1]  used  series  ex¬ 
pressions  in  place  of  (5).  These  are  unsuitable  because  they 
cannot  be  applied  to  all  geometries.  Modem  “math  packs” 
available  on  large  main-frame  computers  can  evaluate  elliptic 
integrals  readily  in  all  situations.  Thus  (S)  was  employed  to 
compute  the  mutual  inductances. 

To  calculate  the  mutual  inductance  between  two  rings  /  and  / 
in  our  model,  both  of  the  rings  are  subdivided  into  n  smaller 
rings.  Each  smaller  ring  is  then  replaced  by  a  filament  in  its 
center.  Following  Chang  [1] ,  the  mutual  inductance  is  then 
given  by 


t  "u-p  <*> 

"  1/  -  l  my- I 

where  Miimj  is  the  mutual  inductance  as  given  by  (5)  of  two 
filaments  lt  and  mf  inside  rings  /  and  /,  respectively. 

To  compute  the  self-inductance  of  a  conducting  ring  /  in  our 
model,  the  ring  is  subdivided  into  n  smaller  rings,  and  these 
smaller  subrings  are  again  replaced  by  filaments.  The  self¬ 
inductance  of  the  conducting  ring  is  then  given  by 


Inn  1  n 

n  1,  -1  mj-l  "  1,-1 


(0) 


(9) 


where  l\^  is  the  self-inductance  of  filament  llt  given  by 
Smythe  [13]  as 


40)  a  n0a,  [In  (16a,/c,)-  1.75]  (10) 


for  a,- »  c,,  where  at  is  the  radius  of  the  filament  and  ct  is  the 
width  of  the  subring  replaced  by  the  filament.1 

The  mutual  inductance  Mt0  between  a  conducting  ring  and 
the  excitation  coil  is  calculated  in  the  same  way  as  other  mu¬ 
tual  inductances.  The  ring  is  divided  into  n  subrings,  the  sub¬ 
rings  are  replaced  by  filaments,  and  the  excitation  coil  of  n0 
turns  is  replaced  by  a  vertical  stack  of  n0  concentric  circular 
filaments,  equally  spaced.  A  multilayered  excitation  coil  is 
modeled  as  several  concentric  vertical  stacks  of  filaments,  ap¬ 
propriately  spaced.  A liQ  is  given  by 


Mto 


» o 


£  Mlimr 

mi  *  I 


The  resistance  of  the  rings  is  calculated  by 


Ri  ~ 


(ID 


(12) 


where  Rt  is  the  resistance  of  the  ith  ring,  p  is  its  resistivity, 
Li  -  2itr(  is  the  mean  path  length  of  current  in  the  ring,  rt  is 
the  mean  radius  of  the  ring,  and  At  is  the  cross-sectional  area 
of  the  ring. 

Current  densities  in  each  ring  are  easily  derivable  from  ma¬ 
trix  [/]  in  (4)  and  from  ring  dimensions. 


'Allhough  this  condition  is  not  satisfied  for  the  innermost  filament 
of  the  innermost  ring,  it  docs  not  lead  to  significant  error  in  the  numeri¬ 
cal  calculation  because,  as  cited  in  the  text,  the  current  will  be  negligi¬ 
ble  in  the  region  where  the  approximation  breaks  down. 


III.  Comparison  to  Finite  Element  Results  for 
the  Case  of  the  Double-Layered  Plate 

Although  we  have  constructed  a  computer  program  to  han¬ 
dle  a  problem  as  general  as  that  depicted  in  Fig.  1,  we  here  re¬ 
strict  the  geometry  to  a  single-layer  20-tum  excitation  coil 
above  a  conducting  plate  of  two  layers.  In  the  first  layer  of 
the  plate,  the  resistivity  p  is  6.6934  X  10" 7  ohm-meters  (ap¬ 
propriate  to  stainless  steel);  in  the  second  layer,  p=  1.5221  X 
10"6  (appropriate  to  certain  titanium  alloys).  The  first  layer 
is  0.127  cm  (50  mils)  thick;  the  second  of  unlimited  thickness. 
The  coil  has  a  length  of  0.3175  cm  (125  mils)  and  diameter 
of  0.3175  cm.  The  coil  carries  0.5  A  of  sinusoidally  varying 
current. 

Presented  here  are  the  results  of  two  numerical  solutions  to 
the  problem -one  by  the  coupled  circuit  approach  and  one  by 
the  finite  element  approach. 

In  the  coupled  circuit  approach,  the  first  layer  of  the  plate  is 
subdivided  into  four  sublayers,  each  0.03175  cm  thick,  and 
the  second  layer  is  subdivided  into  nine  sublayers  of  similar 
thickness.  The  fact  that  the  second  layer  is  given  a  finite 
thickness  has  a  negligible  effect  on  the  results  because  of  the 
skin  effect  and  because  the  plate  is  modeled  as  a  disk  of  suffi¬ 
ciently  large  diameter,  namely,  2.655  cm.  Each  layer  of  the 
plate  is  subdivided  into  12  concentric  subrings  of  equal  radial 
thickness.  In  this  way,  both  a  horizontal  and  vertical  distribu¬ 
tion  of  the  eddy  currents  in  the  plate  can  be  obtained. 

Fig.  2  shows  details  concerning  the  way  in  which  Palanisamy 
[3]  did  the  corresponding  finite  element  calculation.  The  coil 
is  modeled  as  a  region  of  uniform  current  density  4.133  X  103 
A/cms,  of  inner  radius  0.1588  cm  and  outer  radius  0.1664  cm. 
Finite  elements  inside  the  conducting  plate  are  not  uniform  in 
size  but  are  smaller  in  the  regions  of  high  eddy  current  density 
3nd  larger  in  regions  of  low  eddy  current  density,  gradually  in¬ 
creasing  in  size  from  high  to  low  density  regions,  as  seen  in  the 
discretization  depicted  in  Fig.  2.  By  contrast,  for  a  single  sinu¬ 
soidal  excitation  frequency,  the  coupled  circuit  approach  does 
not  require  a  similar  refinement  utilizing  graduated  ring  size  at 
least  for  frequencies  in  the  range  from  5  to  2000  kHz. 

Fig.  3  depicts  the  variation  of  the  eddy  current  density  mag¬ 
nitude  with  depth  at  a  particular  radius  from  the  central  axis, 
namely,  0.1664  cm.  A  characteristic  discontinuity  in  current 
density  versus  depth  is  seen  at  the  boundary  between  the  two 
plate  layers.  Both  finite  element  results  (solid  curves)  and  cou¬ 
pled  circuit  results  (dashed  curves)  are  presented.  Agreement 
in  the  frequency  range  12.5-200  kHz  is  extremely  good  de¬ 
spite  the  small  differences  in  the  way  in  which  the  excitation 
coil  is  modeled. 

The  coupled  circuit  computation  was  performed  on  a  CD C 
CYBER  172  computer  which  is  a  64-bit  machine,  whereas  the 
finite  element  computation  was  performed  on  a  VAX  1 1/780, 
which  is  a  32-bit  machine.  The  coupled  circuit  computation 
involved  initial  construction  of  a  1 56  X  1 56  mutual  induc¬ 
tance  matrix  (requiring  about  7.5  min  of  CPU  time).  Once 
formulated,  this  inductance  matrix  however  did  not  have  to  be 
recomputed  for  each  new  frequency.  The  finite  element  tech- 
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Fig.  2.  Air-core  coil  over  a  conducting  slab  and  the  finite  element  discretization  (courtesy  of  R.  Palanisamy ). 


EDDY  CURRENT  DENSITY  (J)  VS  OEPTH 

Fig.  3.  Current  density  in  plate  a;  function  of  depth  for  frequencies 
from  1 2.3  kHz  to  200  kHz.  The  discontinuity  seen  at  depth  of  0.1 27 
cm  (30  mils)  is  due  to  change  in  medium.  Top  plate  layer  is  stainless 
steel  ar.d  bottom  layer  is  titanium  alloy. 


nique  required  the  storage  of  60  000  element  matrix.  Both 
techniques  required  roughly  one  minute  per  frequency  to  com¬ 
pute  respective  spatial  distributions  of  eddy  currents.  In  terms 
of  computation  times,  the  two  techniques  are  roughly  com¬ 
parable,  with  finite  element  possibly  being  slightly  faster  in  ini¬ 
tial  construction  of  matrices,  but  with  finite  element  also  re¬ 
quiring  slightly  more  storage.  This  latter  however  may  not 
necessarily  be  true  as  perhaps  a  coarser  mesh  could  have  given 
somewhat  similar  results. 

IV.  Comparison  to  Analytical  Results  fox 
the  Case  oi  the  Double-Layered  Plate 
Dodd  el  al.  [7] -[9]  have  found  an  analytical  solution  for 
the  problem  of  a  sinusoidally  varying  current  in  an  excitation 
coil  above  a  double-layered  plate,  where  the  second  layer  is  of 
infinite  vertical  extent  and  both  layers  are  of  infinite  horizon¬ 
tal  extent. 

Two  models  can  be  used  for  the  excitation  coil.  First  the 
coil  can  be  represented  as  a  stack  of  equally  spaced  circular 
filaments  ("delta-function  coils").  Second,  it  can  be  repre¬ 
sented  as  a  cylindrical  sheet  with  length  and  radius  correspond¬ 
ing  to  that  of  the  excitation  coil  and  with  current  per  unit 
length  equal  to  ft0/0/(/„0  -  /,),  where  /,  and  l»o  are  distances 
from  the  plate  to  the  bottom  and  top  of  the  excitation  coil, 
respectively,  and  where  the  coil  has  n0  turns. 
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In  the  case  of  the  first  model,  the  current  densities  and 
7,(/)  in  plate  tayers  1  and  2.  respectively,  are  given  in  phasor 
notation  by  {7] -[9J 


(f) 


('.*) 


f  Jl 


1  t  -a0l  i 


)mo,ji0/0'o  J  y,(nr0)y,(ar)  - - e 


l-e 


X  a 


4'V.z) 


(g,  ]  da 


03a) 


Jr-  j  _  e-"o «oAt 

!  ■M«',©)-M®') - e 

o  1  -  e 


J* i  +®i)0  «:« 

_ *P|g _ £ _ 

(«o  -  0, X0.  fc)  +  (*o  ♦  0i W . ♦  ft)  e10'" 

03b) 

where  o(  and  o,  are  the  respective  conductivities  of  the  two 
plate  layers,  r0  is  the  radius  of  the  excitation  coil,./|(ar0)  is  a 
Bessel  function  of  order  one,  A/  3  (f„0  -  /,  )/n0,D  is  the  thick¬ 
ness  of  the  first  plate  layer,  and 

a,  3  V<*2  -  u'njej  +jwn/Oj  ,  (14a) 

ft=(Po/P/)a/.  (Hb) 


with  j  3  0, 1 , 2  and  o0  3  0.  The  numerical  results  presented 
earlier  have  been  calculated  under  the  assumption  that  p0  3 
H{  3  Mj  and  e0  =  Ci  3  e2l  so  that  =  otj.  For  the  case  of  the 
second  model,  the  expression 


I  .  e-»o<»oA' 


1  -  e 


-ao  A/ 


e'ao'i 


(IS) 


is  replaced  by 


The  first  analytical  model  shall  be  called  the  “g-model,”  since 
a  geometric  series  is  used  in  the  evaluation.  The  second  ana¬ 
lytical  model  shall  be  called  the  “/-model,”  since  an  integral 
over  the  length  of  the  coil  is  used.  Two  models  are  introduced 
because  in  the  coupled  circuit  analysis,  the  coil  was  modeled  as 
a  stack  of  circular  filaments  discretely  spaced  (which  corre¬ 
sponds  to  the  g-model),  whereas  in  the  finite  element  case,  the 
coil  was  treated  as  a  narrow  cylinder  of  constant  current  den¬ 
sity  (which  corresponds  roughly  to  the  /  model). 

Fig.  4  presents  the  results  of  both  the  g-model  and  the  /• 
model  for  the  geometry  described  in  Section  111.  The  coupled 
circuit  results  are  superimposed  on  the  analytic  curves.  It  is 
seen  that  agreement  is  extremely  good.  The  slight  differences 
between  g-model  results  and  /-model  results  are  about  the  same 
as  the  differences  between  the  coupled  circuit  results  and  the 
finite  element  results  noted  earlier,  suggesting  that  the  way  in 
which  the  coil  is  modeled  chiefly  accounts  for  differences  be¬ 
tween  results  from  the  two  numerical  approaches. 


DEPTH  Z  (cm) 

Fig.  4.  Current  density  in  two-layered  plate  as  function  of  depth.  Solid 
curves  represent  analytical  results,  with  label  g  denoting  “y-modcl” 
and  label  /  denoting  “/-model,"  as  explained  in  text.  Solid  circles 
show  the  results  from  coupled  circuit  computation. 

Figs.  5  and  6  focus  on  a  particularly  interesting  effect,  namely, 
that  the  current  density  at  the  boundary  between  the  two 
plate  layers  exhibits  a  broad  maximum  in  its  frequency  depen¬ 
dence  for  the  range  300-500  kHz.  The  coupled  circuit  ap¬ 
proach  and  the  analytic  approach  both  exhibit  this  same  ef¬ 
fect,  whereas  the  finite  element  calculation  was  not  tried  for 
this  frequency  range. 

Fig.  7  shows  a  comparison  between  analytic  and  coupled  cir¬ 
cuit  approaches  for  a  frequency  range  in  order  of  magnitude 
greater  than  that  seen  in  Fig.  3.  As  can  be  seen,  the  coupled 
circuit  approach  gives  fairly  good  predictions  without  an  ad¬ 
justment  in  ring  subdivisions. 

V.  Cknkralization  to  Nonsinusoidal  Excitation 

The  coupled  circuit  approach  can  be  extended  to  nonsinu¬ 
soidal  excitation,  but  at  the  expense  of  considerable  computa¬ 
tion  time.  The  first  step  is  to  obtain  a  fast  Fourier  transform 
of  the  excitation  signal  /<>(/).  Thus,  for  example,  values  /0(w*) 
are  obtained  for  the  discrete  set  of  frequencies  {w*}.  The  re- 
sponse  lt(uk)  of  the  ith  conducting  ring  to  sinusoidal  excita¬ 
tion  current  /0<',w*r  is  then  obtained  using  coupled  circuit 
solution  (4).  A  fast  Fourier  transform  is  then  taken  of  the 
product  set  {/,,( co*)  //(w*)}  to  obtain  the  time-developing 
current  /,(/)  in  ring  i.  A  similar  numerical  transform  proce¬ 
dure  for  pulsed  excitations  has  been  reported  recently  by 
Allen  { 1 4 j . 
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DEPTH  Z  (cm) 

Fig.  S.  Current  density  versus  depth  for  higher  frequencies  in  range  from  100-500  kHz.  Solid  curves  denote  /-model  ana¬ 
lytical  results;  dashed  curves,  g-model  analytical  results.  Coupled  circuit  results  are  denoted  by  •,  100  kHz;  o,  300  kHz; 
and  o,  500  kHz.  Crasser  denote  300  kHz  p-modeTresults  at  depths  where  300  kHz  and  500  kHz  curves  intersect. 


FREQUENCY  (KHZ) 

Fig.  6.  Current  density  versus  frequency  at  depth  of  0.127  cm,  corresponding  to  boundary  between  layer  I  (stainless  steel) 
and  layer  2  (titanium  alloy).  Analytical  results  are  given  by  solid  curves.  Top  two  curves  represent  currents  on  stainless 
steel  tide  of  the  boundary;  bottom  curves,  the  titanium  side. 
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1  -  1000  KHl 

2  -  2000  KHl 


3  400 


OEPTH  Z  (cm) 

Fig.  7.  Current  density  versus  depth  for  frequencies  1000  kHz  (•)  said 
2000  kHz  (o).  As  before,  solid  curves  represent  analytical  results 
(labeled  f  and  /),  an d  points  represent  coupled  circuit  results. 


We  have  tried  this  approach  on  the  problem  of  an  excitation 
coil  surrounded  by  a  conduction  shield.  This  problem  has  rele¬ 
vance  to  probe  design  in  the  pulsed  eddy  current  nondestruc¬ 
tive  evaluation  technique  [IS].  An  analytical  solution  is  now 
available  for  this  situation  [IS]. 

For  this  case,  uniform  discretization  into  conducting  rings 
does  not  produce  results  that  completely  represent  the  correct 
time  development  of  the  current  response.  Much  better  results 
are  obtainable  if  rings  are  narrow  and  thin  in  regions  of  high 
current  density  and  broad  and  thick  in  regions  of  low  current 
density.  Thus  for  nonsinusoidal  excitation,  nonuniform  dis¬ 
cretization  is  necessary. 
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VI.  Conclusion 

An  alternative  numerical  approach  to  the  problem  of  com¬ 
puting  induced  eddy  current  densities  in  axisymmetric  geome¬ 
tries  has  been  discussed.  This  “coupled  circuit  approach"  re¬ 
produces  analytical  solutions  extremely  well  in  the  case  of 
single-frequency  excitation  using  a  relatively  simple  discretiza¬ 
tion  procedure.  Computation  times  and  storage  requirements 
are  comparable  to  the  finite  element  approach.  The  coupled 
circuit  approach  can  also  be  used  for  nonsinusoidal  excitation 
using  a  fast  Fourier  transform  procedure. 
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APPENDIX  C 


PROBE  DESIGN  STUDIES 


Transmitter  Impulse  Response  Calculations 

The  transmitter  impulse  response  function  E/I  is  the  transient  electric 
field  induced  in  the  specimen  by  a  current  pulse  of  zero  width  in  the 
transmitter  coil.  The  purpose  of  the  calculations  reported  here  was  to 
determine  the  effects  of  transmitter  geometry  on  induced  fields,  and  thus 
to  provide  guidance  in  designing  the  transmitter  to  maximize  the  induced 
field  at  typical  flaw  locations. 

Figure  43  shows  the  geometry  used  for  all  of  the  calculations  in  this 
Appendix.  In  some  of  the  calculations,  such  as  those  presented  in  the 
next  four  figures,  the  effect  of  the  flaw  on  the  receiver  response  is  not 
involved.  Thus,  the  only  parameters  of  concern  at  present  are  those  per¬ 
taining  to  the  transmitter  geometry  and  position  (p,  Z)  within  the  speci¬ 
men.  Values  chosen  for  the  initial  study  were  the  following: 


Coil  length  (LT) 

Core  radius  (CT) 
Winding  thickness  (W-j) 
Shield  radius  (ST) 
Lift  off  (ZT) 


-  3.2  mm  (0.125”) 

-  1.02  mm  (0.040") 

-  1.01  mm  (0.040”) 

-  2.29  mm  (0.090”) 
*  0.42  mm  (0.016”) 


Figures  44  through  48  show  the  effects  of  variations  In  transmitter  geom¬ 
etry  on  the  induced  field.  Ferrite  core  diameter  is  an  important  param¬ 
eter,  as  illustrated  in  Figure  44,  which  shows  how  the  Induced  field  is 
broadened  and  increases  in  intensity  as  core  size  Increases.  Variations 
in  coll-to-shield  distance,  on  the  other  hand,  produce  only  minor  changes 
in  the  field,  as  shown  in  Figure  45.  The  next  illustration.  Figure  46, 
illustrates  the  expected  lift-off  effect;  and  the  final  figure.  Figure  47, 
shows  how  the  shield  tends  to  focus  the  Induced  field,  again  as  expected. 


Flaw  Impulse  Response  Function 

The  flaw  impulse  response  function  is  proportional  to  the  change  in  volt¬ 
age  in  the  receiver  coil  caused  by  a  flaw,*  with  the  source  current  being 
a  pulse  of  zero  width  in  the  transmitter  coil.  In  the  next  two  figures, 
flaw  response  calculations  are  presented  as  a  function  of  flaw  position, 
receiver  winding  configuration,  and  signal  arrival  time. 

The  detailed  model  treats  a  differential  receiver  with  alternating  posi¬ 
tive  and  negative  windings,  as  illustrated  in  Figure  43,  where  only  one 
pair  of  windings  is  shown.  To  study  the  effect  of  winding  configuration, 
calculations  were  performed  as  a  function  of  the  number  of  ±  winding  pairs 
with  the  overall  thickness  of  the  winding  held  constant.  Figure  48  shows 
a  typical  comparison,  in  this  case  between  a  coil  with  one  ±  layer 

*In  these  calculations,  as  in  Burrows'  (9)  model,  the  flaw  is  assumed  to 
be  small  enough  to  be  treated  as  a  point  scatterer.  Signal  amplitudes  are 
normalized  by  dividing  by  the  flaw  volume. 


Figure  43.  Geometry  Used  in  Probe  Design  Studies 
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(NLAYERS  *  1)  and  a  coil  with  5  such  layers  (NLAYERS  *  5).  As  can  be  seen 
from  these  data,  the  principal  effect  of  increasing  the  number  of  layers 
is  to  decrease  the  signal  voltage  (the  NLAYERS  -  5  curve  has  been  multi¬ 
plied  by  5)  with  no  improvement  in  flaw  resolution.  Another  point  illus¬ 
trated  here  is  that  flaw  response  is  greatest  when  the  flaw  is  directly 
under  the  transmitter  winding  on  the  side  nearest  the  detector  (ps  1.0  mm) 
in  agreement  with  experimental  observations.  Similar  calculations  at  dif¬ 
ferent  flaw  depths  show  that  the  shape  of  the  flaw  response  curve,  as 
plotted  in  Figure  48,  does  not  change  significantly  with  flaw  depth. 

All  calculations  of  the  time-dependent  signal  produced  in  the  receiver 
coil  show  certain  common  features,  which  are  depicted  in  Figure  49.  The 
impulse  response  signal  is  characterized  by  a  very  sharp  positive  peak 
followed  by  a  series  of  secondary  peaks  of  decreasing  amplitude,  increas¬ 
ing  width,  and  alternating  sign.  Peak  amplitudes  and  widths  depend  pri¬ 
marily  on  flaw  depth  and  are  not  strongly  influenced  by  p,  the  radial  dis¬ 
tance  from  the  transmitter  to  the  flaw.  This  behavior  is  illustrated  by 
data  presented  in  Figure  46  of  the  text  for  a  short  Gaussian  excitation 
pulse.  Peak  arrival  times  showing  the  depth  dependence  are  illustrated  in 
Figure  47. 


Figure  49.  Schematic  Illustration 
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APPENDIX  D 


VOLTAGE  TO  CURRENT  CONVERTER 


The  function  of  this  circuit  is  to  produce  an  output  current  whose  ampli¬ 
tude  is  proportional  to  the  amplitude  of  the  input  voltage.  Figure  50  is 
a  schematic  diagram  of  the  final  version  adapted  for  the  pulse  shaping 
experiments.  Following  is  a  brief  description  of  the  circuit. 

The  input  voltage  is  amplified  by  amplifer  AR1  and  subsequently  inpressed 
on  the  base  of  transistor  Ql.  In  the  absence  of  a  feedback  current  through 
resistor  R13,  the  gain  of  AR1  is  about  10;  however,  when  an  appropriate 
load  is  connected  to  the  output  of  the  circuit,  a  feedback  current  is 
generated  that  reduces  the  effective  gain  to  about  5. 

Transistor  Ql  acts  to  convert  the  positive  going  voltage  at  its  base  into 
a  negative  going  current  in  its  collector  circuit.  The  resistors  connected 
in  the  emitter  circuit  perform  a  degenerative  function  that  acts  to  make 
the  conversion  linear.  The  capacitor  connector  to  the  emitter  enhances 
the  high  frequency  response  of  transistor.  The  collector  current  passes 
through  resistor  R7  to  produce  a  negative  going  voltage  at  the  input  of 
buffer  amplifier  AR2.  AR2  is  a  unity  gain,  non-inverting  amplifier  with 
considerable  drive  capability.  Its  output  drives  the  high  susceptance 
input  gate  of  the  HEXFET  transistor  Q2.  The  maximum  drive  voltage 
required  at  this  point  is  about  10  volts  (negative). 

The  HEXFET  produces  a  current  that  is  a  function  of  the  drive  voltage  at 
its  gate.  This  current  passes  through  the  device  connected  between  the 
output  terminals  and  through  the  2  ohm  resistor  R14.  The  output  current 
of  the  HEXFET  is  not  a  linear  function  of  its  drive  voltage.  To  correct 
for  this  non-linearity  and,  others  that  might  be  in  the  signal  circuit, 
the  voltage  developed  across  R14  is  used  to  produce  a  feedback  current 
through  R13.  This  current  acts  to  modify  the  output  of  AR1  in  a  manner 
that  will  compensate  for  the  circuit  non-linearities  and  make  the  ampli¬ 
tude  of  the  output  current  a  proportional  replica  of  the  amplitude  of  the 
input  voltage. 

The  circuit  parameters  have  been  chosen  such  that  a  two  volt  input  pulse 
produces  an  output  current  pulse  of  approximately  2  amperes.  The  loop 
gain  of  the  feedback  circuit  has  been  adjusted  to  give  the  nearest  approach 
to  a  linear  relationship  between  the  input  voltage  and  the  output  current 
that  can  be  achieved  without  danger  of  oscillation  or  significant  transient 
ringing.  Better  results  might  be  achieved  if  the  circuit  is  redesigned  to 
use  higher  speed  transistors. 

NOTES: 

(1)  The  components  connected  in  the  source  circuit  of  Q3  are  carry-overs 
from  the  original  circuit  design  and  are  of  insignificant  importance 
in  its  current  application. 


(2)  The  transistor  Q4  acts  to  provide  a  small  bias  current  that  over¬ 
comes  and  leakage  current  that  might  be  present  and  aasures  that  the 
voltage  at  the  output  terminals  is  zero  under  no  drive  conditions* 

(3)  The  transistors  Q3  and  Q6  act  with  associated  circuitry  to  limit  the 
high  voltage  current  demand  to  approximately  20  mA. 

(4)  Because  Q2  acts  as  a  current  source,  the  value  of  the  high  voltage 
has  no  significant  effect  on  the  operation  of  the  circuit. 


Figure  50.  Schematic  of  Voltage  tp  Current  Converter 
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